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We report results from a combined analysis of solar neutrino data from all phases of the Sudbury 
Neutrino Observatory. By exploiting particle identification information obtained from the propor- 
tional counters installed during the third phase, this analysis improved background rejection in that 
phase of the experiment. The combined analysis resulted in a total flux of active neutrino flavors 
from *B decays in the Sun of (5.25 ± 0.16(stat.)^Q' ^^(syst.)) x 10*^ cm~^s~^. A two-flavor neutrino 
oscillation analysis yielded Amji = (5. 4) x lO"'^ eV^ and tan^ 6^12 = 0.427to {jag- A three-flavor 
neutrino oscillation analysis combining this result with results of all other solar neutrino experiments 
and the KamLAND experiment yielded Amii = (7.4llJ5 ?9) x 10"^ eV^ tan^ 612 = 0.446l2;^|^, and 
sin'^ 613 = {2.5t-i's) X lO"'^. This implied an upper bound of sin^ 613 < 0.053 at the 95% confldence 
level (C.L.). 



I. INTRODUCTION decays in the Sun, so-called *B neutrinos, and to study 

neutrino oscillations, as proposed by Serb Chen [T]. As a 
The Sudbury Neutrino Observatory (SNO) was de- result of measurements with the SNO detector and other 
signed to measure the flux of neutrinos produced by ^B 
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experiments, it is now well-established that neutrinos are 
massive and that the weak eigenstates {vei v^^, Vr) are 
mixtures of the mass eigenstates (i^i, i^2: '^a)- The prob- 
ability of detecting a neutrino in the same weak eigen- 
state in which it was created depends on the energy and 
propagation distance of the neutrino, the effects of mat- 
ter [H [3], the neutrino mixing angles (^12, 6*23, ^13), a 
phase (S) which can lead to charge-parity violation, and 
the differences between the squares of the neutrino mass 
eigenvalues {Am^i, Ato|2, Am|]^) [4 



32, ^'"-31 

The SNO detector observed ®B neutrinos via three dif- 
ferent reactions. By measuring the rate of neutral current 
(NC) reactions, 

+ d ^ p + n + (1) 

which is equally sensitive to all three active neutrino 
flavors, the SNO experiment determined the total 
neutrino flux, $Bj independently of any specific active 
neutrino flavor oscillation hypothesis [T]. The predicted 
flux from solar model calculations [B] is (5.88 ± 0.65) x 
10^ cm-2s-i, BPS09(GS), or (4.85±0.58) x 10^ cm-^g-i, 
BPS09(AGSS09), using a recent measurement of the 
heavy-element abundance at the Sun's surface. Previ- 
ous analyses of SNO data [3 [5j measured $b more pre- 
cisely than the solar model predictions. A more precise 
measurement of <i>B would better constrain these solar 
models, but may not necessarily determine which metal- 
licity is correct due to the large uncertainties at present 
on both predictions. 

By measuring the rate of charged current (CC) reac- 
tions, 

Ve+d^p+p + er, (2) 

which is only sensitive to z^gS, and comparing this to the 
NC reaction rate, it was possible to determine the neu- 
trino survival probability as a function of energy. This 
can then constrain the neutrino oscillation parameters 
independently of any specific prediction of $3 ■ 

The SNO experiment also measured the rate of elastic 
scattering (ES) reactions, 

i^x + ^ i^x + e~ , (3) 

which is sensitive to all neutrino fiavors, but the cross- 
section for VeS is approximately six times larger than that 
for the other fiavors. 

We present in this article a final combined analysis of 
all solar neutrino data from the SNO experiment. Sec- 
tion [n] gives an overview of the detector. In Section III 
we describe the method used to combine all the data 
in a fit which determines <1>b and a parameterized form 
of the i^e survival probability. Section |IV| describes a 
new particle identification technique that allowed us to 
significantly suppress the background events in the pro- 
portional counters used in the third phase of the SNO 
experiment. Section |V] presents the results of the new 
analysis of data from Phase III, and the combined analy- 
sis of data from all phases. The results of this combined 
analysis are interpreted in the context of neutrino oscil- 
lations in Section IVll 




FIG. 1. Schematic diagram of the SNO detector. We used 
a coordinate system with the center of the detector as the 
origin, and z direction as vertically upward. 



II. THE SNO DETECTOR 

The SNO detector [3] , shown schematically in Figure 
consisted of an inner volume containing 10® kg of 99.92 
isotopically pure heavy water (^Il20, hereafter referred 
to as D2O) within a 12 m diameter transparent acrylic 
vessel (AV). Over 7x 10® kg of II2O between the rock and 
the AV shielded the D2O from external radioactive back- 
grounds. An array of 9456 inward-facing 20 cm Hama- 
matsu R1408 photomultiplier tubes (PMTs), installed 
on an 17.8 m diameter stainless steel geodesic structure 
(PSUP), detected Cherenkov radiation produced in both 
the D2O and H2O. A non-imaging light concentrator [TU] 
mounted on each PMT increased the effective photocath- 
ode coverage to nearly 55% of 47r. The PMT thresholds 
were set to 1/4 of the charge from a single photoelectron. 
The inner 1.7 X 10® kg of H2O between the AV and the 
PSUP shielded the D2O against radioactive backgrounds 
from the PSUP and PMTs. Extensive purification sys- 
tems removed radioactive isotopes from both the D2O 
and the H2O [II]. 

The detector was located in Vale's Creighton mine 
(46°28'30"N latitude, 81° 12'04" W longitude) near Sud- 
bury, Ontario, Canada, with the center of the detector 
at a depth of 2092 m (5890±94 meters water equivalent). 
At this depth, the rate of cosmic-ray muons entering the 
detector was approximately three per hour. Ninety-one 
outward-facing PMTs attached to the PSUP detected 
cosmic-ray muons. An offiine veto based on informa- 
tion from these PMTs significantly reduced cosmogenic 
backgrounds. 

The recoil electrons from both the ES and CC reac- 
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tions were detected directly through their production of 
Cherenkov hght. The total amount of light detected by 
the PMT array was correlated with the energy of the 
interacting neutrino. 

The SNO detector operated in three distinct phases 
distinguished by how the neutrons from the NC inter- 
actions were detected. In Phase I, the detected neu- 
trons captured on deuterons in the D2O releasing a sin- 
gle 6.25 MeV 7-ray, and it was the Cherenkov light of 
secondary Compton electrons or e~^e~ pairs that was 
detected. In Phase II, 2 x 10"^ kg of NaCl were added 
to the D2O, and the neutrons captured predominantly 
on ^^Cl nuclei, which have a much larger neutron cap- 
ture cross-section than deuterium nuclei, resulting in a 
higher neutron detection efficiency. Capture on chlorine 
also released more energy (8.6 MeV) and yielded multi- 
ple 7-rays, which aided in identifying neutron events. In 
Phase III, an array of proportional counters (the Neutral 
Current Detection, or NCD, array) was deployed in the 

D2O [H]. 

The proportional counters were constructed of approx- 
imately 2 m long high purity nickel tubes welded together 
to form longer "strings" . Neutrons were detected via the 
reaction 



^He + n ^ -^R + p, 



(4) 



where the triton and proton had a total kinetic energy of 
0.76 MeV, and travelled in opposite directions. The NCD 
array consisted of 36 strings filled with '^He, and an ad- 
ditional 4 strings filled with ^He that were insensitive to 
the neutron signals and were used to study backgrounds. 
Energetic charged particles within the proportional coun- 
ters produced ionization electrons, and the induced volt- 
age caused by these electrons was recorded as a function 
of time, referred to as a waveform. To increase the dy- 
namic range, the waveform was logarithmically amplified 
before being digitized [T^ . 



the Sun, cos^©. Different algorithms to calculate both 
Toff and p were used for the first two phases and the third 
phase. References |71[T31[I1] contain detailed descriptions 
of how these variables were calculated. The energy de- 
posited in the gas of a proportional counter, i?NCD, was 
calculated for each event recorded with the NCD array, 
and the correlated waveform was determined [14 . 

Although there were multiple sets of data in this fit, the 
result was a single $b and energy-dependent survival 
probability as described in Section [ill. 1[ We summarize 
the event selection and backgrounds in Section [ill. 2 [ Sec- 
tions [ITL3| and |nL4| respectively, describe the PDFs and 
efficiencies. The method for combining the multiple sets 
of data in a single analysis is presented in Section [III.5[ 
Finally, Section p^II.6[ outlines the alternative analyses to 
verify the combined analysis. 



III.l. Parameterization of the neutrino signal 

We fitted the neutrino signal based on an average $b 
for day and night, a survival probability as a function 
of neutrino energy, E^, during the day, P^^{E^), and an 
asymmetry between the day and night survival probabil- 
ities, Aee{E^), defined by 



Aee{E,) = 2 



P-,{E,) + pi,{E,y 



(5) 



where P^^(E^) was the i^e survival probability during the 
night. This was the same parameterization as we used in 
our previous analysis of data from Phases I and II [7]. 

Monte Carlo simulations assuming the standard solar 
model and no neutrino oscillations were used to deter- 
mine the event variables for ®B neutrino interactions in 
the detector. These simulations were then scaled by the 
factors given in Table [l[ 



III. COMBINED ANALYSIS 

In this article we present an analysis that combines 
data from all three phases of the SNO experiment. The 
combination accounts for any correlations in the system- 
atic uncertainties between phases. The data were split 
into day and night sets in order to search for matter ef- 
fects as the neutrinos propagated through the Earth. 

The general form of the analysis was a fit to Monte 
Carlo-derived probability density functions (PDFs) for 
each of the possible signal and background types. As 
with previous analyses of SNO data, the following four 
variables were calculated for each event recorded with the 
PMT array: the effective electron kinetic energy, T^s, 
reconstructed under the hypothesis that the light was 
caused by a single electron; the cube of the radial posi- 
tion, r, divided by 600cm, p — {^[™]/60o)^; the isotropy 
of the detected light, /3i4; and the angle of the recon- 
structed electron propagation relative to the direction of 



TABLE I. **B neutrino interactions scaling factors. 



Interaction Day/Night 


Scaling 


factor 


CC,ESc 


Day 








Day 


*b[1 - 




CC,ESe 


Night 






ES^x 


Nigiit 


*b[1 - 


PfAE^)] 


NC 


Day+Niglit 


*B 





Unlike our earlier publications [8j [131 [E] : this param- 
eterization included a constraint on the rate of ES inter- 
actions relative to the rate of CC interactions based on 
their relative cross-sections. It also had the advantage 
that the fitted parameters ($b, PteiEy)^ and A^eiEu)) 
were all directly related to the scientific questions under 
investigation. Moreover, it disentangled the detector re- 
sponse from the fit result as P^^{E^) and Aee{E^) were 
functions of E^, as opposed to Toff. 

Appendix [^ explains how this parameterization can 
be used to describe sterile neutrino models that do not 



4 



predict any day/night asymmetry in the sterile neutrino 
flux and do not predict any distortion in the sterile 
spectrum [12]. Reference [T7] presents a very general 
sterile neutrino analysis that includes day versus night 
asymmetries. 

Due to the broad Toff resolution of the detector, 
Pf^{E^) was not sensitive to sharp distortions and was 
parameterized by 



Pi{E,) = Co + ci(^,[MeV] - 10) 
+C2(S,[MeV] - 10)2, 



(6) 



where cq , ci , and C2 were parameters defining the Ve sur- 
vival probability. Simulations showed that the fit was not 
sensitive to higher order terms in the polynomial. Ex- 
panding the function around 10 MeV, which corresponds 
approximately to the peak in the detectable *B neutrino 
En spectrum, reduced correlations between cq, ci, and 
C2. For the same reasons, Aee{Ei,) was parameterized by 



AeeiE,) =ao + ai(S,[MeV] - 10), 



(7) 



where ap, and ai were parameters defining the relative 
difference between the night and day survival prob- 
ability. By disallowing sharp changes in the neutrino 
signal that can mimic the background signal at low en- 
ergies, these parameterizations reduced the covariances 
between the neutrino interaction and background rates. 

To correctly handle ES events, we simulated v^s with 
the same E^, spectrum as VeS, such that scaling factors 
for these interactions in Table[l]were satisfied. In our pre- 
vious analysis [7] we approximated the and Vr cross- 
section by 0.156 times the Ve cross-section, and then in- 
cluded an additional systematic uncertainty to account 
for the fact that the ratio of the I'g to j/^ ES cross-section 
is not constant as a function of E,,. 



III. 2. Event selection and backgrounds 

Table [n] summarizes the data periods used in this anal- 
ysis. We used the same periods of data as our most recent 
analyses of data from these phases [71 18] . 



TABLE 11. Dates for the data in the different phases used in 
this analysis. 



Phase 


Start date 


End date 


Total time 


[days] 








Day 


Night 


I 


November 1999 


May 2001 


119.9 


157.4 


II 


July 2001 


August 2003 


176.5 


214.9 


III 


November 2004 


November 2006 


176.6 


208.6 



Event cuts to select good candidates were identical 
to those in the previous analyses of these data [3 [S]. 
The following cuts on the event variables were applied: 
p < (550[cm]/600[cm])^ ^ 0.77025, -0.12 < /3i4 < 0.95, 
3.5 MeV < Teff < 20.0 MeV for Phases I and II, and 
6.0 MeV < Toff < 20.0 MeV for Phase III. After these 



cuts the data consisted of events from ES, CC, and NC 
interactions of neutrinos, and a number of different 
background sources. 

Radioactive decays produced two main background 
types: "electron-like" events, which resulted from (3- 
particles or 7-rays with a total energy above our 
Toff threshold, and neutrons produced by the photo- 
disintegration of deuterons by 7-rays with energies above 
2.2 MeV. During Phase III, only the neutron events were 
observed from radioactive background decays, due to the 
higher Tcff threshold for that phase. 

The radioactive decays of ^^^Bi and ^os^pj .^vithin the 
regions of the detector filled with D2O and II2O were 
major sources of background events. ^^''Bi is part of the 
2^^U decay chain, but it was most likely not in equilib- 
rium with the early part of the decay chain. The most 
likely source of ^^''Bi was from ^^^Rn entering the D2O 
and H2O from mine air. "^^^Tl is part of the ^■^^Th de- 
cay chain. These sources of radiation produced both 
electron-like events and photo-disintegration neutrons. 
Ex- situ measurements |18|,I19| of background levels in the 
D2O and II2O provided constraints on the rate of these 



decays, as given in Tables XVIII and XIX| of AppendixfB] 
Background sources originating from the AV included 
decays of ^os^j^ within the acrylic, which produced both 
electron-like events and photo-disintegration neutrons. 
In addition, radon progeny that accumulated on the sur- 
face of the AV during construction could create neutrons 
through (a,n) reactions on isotopes of carbon and oxy- 
gen within the AV. Near the T^ff threshold in Phases I 
and II the majority of background events originated from 
radioactive decays within the PMTs. 

Due to the dissolved NaCl in the D2O during Phase 
II, calibration sources that produced neutrons, and other 
sources of neutrons, led to the creation of ^^Na via 
neutron capture on ^'^Na. ^^Na decays with a half-life 
of approximately 15 hours, producing a low energy (3 
particle and two 7-rays. One of these 7-rays has an 
energy of 2.75 MeV, which could photo-disintegrate a 
deuteron. This resulted in additional electron-like events 
and photo-disintegration neutrons during Phase II. Peri- 
ods after calibrations were removed from the data, but 
there were remaining backgrounds. 

During Phase III there were additional photo- 
disintegration neutron backgrounds due to radioactivity 
in the nickel and cables of the NCD array, as well as two 
"hotspots" on the strings referred to as K2 and K5. The 
estimated number of these background events, given in 
Table |XIX| of Appendix [B] were based on previous anal- 
yses of these data 14, except for backgrounds from the 
K5 hotspot, which was based on a recent reanalysis |20j . 
The previously estimated number of neutrons observed 
in the NCD array due to the K5 hotspot was 31.6 ± 3.7, 
which assumed there was Th and a small amount of U in 
the hotspot based on both ex-situ and in-situ measure- 
ments. Based on measurements performed on the strings 
after they were removed from the D2O, it was determined 
that the radioactivity was likely on the surface and most 
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likely pure Th with very little U. This resulted in a new 
estimate of 45.5lg'4 neutron background events observed 
in the NCD array from this hotspot. 

Aside from the radioactive decay backgrounds, there 
were additional backgrounds to the neutrino measure- 
ment due to j/gS produced by the following reaction, 

^He + p^ ''Hc + e+ + i/e, (8) 

in the Sun, so-called hep neutrinos. These have a maxi- 
mum energy of 18.8 MeV, which is slightly above the ®B 
neutrino maximum energy of 15 MeV, and the standard 
solar model (SSM) prediction for their flux is approxi- 
mately one thousand times smaller than $b pP. Esti- 
mates of this background were based on the SSM pre- 
diction including the effects of neutrino oscillations ob- 
tained from previous analyses [7] . There were instrumen- 
tal backgrounds that reconstructed near the AV. Above 
Toff ~ 6 MeV these events formed a distinct peak at low 
values of /3i4, so they were easily removed by the cuts on 
/3i4 and p. At lower Teff, position reconstruction uncer- 
tainties increase, and the /3i4 distribution of these "AV 
instrumental background" events broadens, resulting in 
incomplete removal by these cuts. This background was 
negligible in Phase III due to the higher Tog threshold 
used for the analysis of data from that phase. Finally, 
there were also background events due to neutrinos pro- 
duced by particle decays in the atmosphere. The esti- 
mated numbers for these background events, given in 
Tables |XVIII| and pClX] of Appendix |B] were based on 
previous analyses of these data \7l ll4j. 

III.3. PDFs 

For Phases I and II the event variables T^s, p, Pi4, 
and cos 00 were used to construct 4-dimensional PDFs. 
For Phase III the reduced number of NC events observed 
with the PMT array made the /3i4 event variable unnec- 
essary, so the PDFs were 3-dimensional in the remaining 
three event variables. Monte Carlo simulations were used 
to derive the PDFs for all signal and background classes 
observed with the PMT array except for backgrounds 
originating from radioactivity in the PMTs, which was 
described by an analytical function. Compared to the 
previous analysis of data from Phases I and II [7], we 
increased the number of Monte Carlo events for the CC 
and ES interactions by a factor of four, and for NC inter- 
actions and some background types by a factor of two. 

The Monte Carlo simulation was verified using a vari- 
ety of calibration sources. Based on these comparisons a 
number of systematic uncertainties were defined to repre- 
sent possible variations in the event variables relative to 
the calibrations. In general these included differences in 
the offset, scale, and resolution for each of the event vari- 
ables. Appendix [B] gives the complete list of systematic 
uncertainties associated with the PDFs. Except where 
specified these uncertainties were the same as those used 
in the most recent analyses of these data [71 [T3] . 



Extensive calibrations using a -'^^N 7-ray source |22) . 
which produced electrons with kinetic energies of approx- 
imately 5.05 MeV from Compton scattering, allowed us 
to calibrate Tog. In Phase I the linearity of T^s with re- 
spect to electron kinetic energy was tested using a proton- 
triton fusion 7-ray source [23], which produced electrons 
with kinetic energies up to approximately 19.0 MeV from 
Compton scattering and pair-production. Based on these 
sources, we parameterized the reconstructed electron ki- 
netic energy including a possible non-linearity by 

/ . reff[MeV]-5.05 \ 
Teff - Teff (^1 + Co ^9 Q_5 05 j, (9) 

where Cq represents the level of non-linearity. The linear- 
ity in all phases was tested using the following two elec- 
tron sources: a ®Li calibration source [21| that produced 
electrons with a continuous distribution up to approxi- 
mately 13 MeV; and electrons with a continuous energy 
distribution up to approximately 50 MeV produced by 
the decay of muons that stopped within the AV. These 
studies revealed non-linearities consistent with zero. We 
assumed any non-linearities below our level of sensitivity 
were correlated between all three phases, and we used a 
value of cf = ± 0.0069, which was equal to the value 
used in the previous analysis of data from Phases I and 
II 7J. 

During Phase III, ^^NaCl brine was injected into the 
D2O on two separate occasions [5S]. The brine was thor- 
oughly mixed in the D2O and provided a uniformly dis- 
tributed source of 7-rays, allowing us to study possible 
Toff variations in regions that were previously not sam- 
pled due to the restricted movement of the -'^^N source. 
The observed variation in the event rate of ^"^Na de- 
cays within the fiducial volume of solar neutrino analysis 
was consistent with what was allowed by the Toff cali- 
bration parameters determined with the ^^N source at 
Toff > 6 MeV. 

While the intrinsic rate of radioactive backgrounds 
from solid bulk materials such as the acrylic vessel or 
PMT array were not expected to vary over the course of 
the experiment, variations in detector response make the 
detected rates vary over time, and because of differences 
in the livetime fractions between day and night, these 
variations were aliased into day/night differences. PDFs 
derived from Monte Carlo naturally include day/night 
detector asymmetries because the detector simulation 
tracks changes to the detector response. Our previous 
analysis derived the analytical PDFs for radioactivity 
originating from the PMTs using a bifurcated analysis 
of real data with the day and night data combined [7], 
which did not account for possible day/night asymme- 
tries. 

To accommodate such asymmetries in the present anal- 
ysis we allowed different observed background rates be- 
tween day and night, and we repeated the bifurcated 
analysis with the data separated into day and night sets. 
Similarly to the previous analysis [71 the PDF was pa- 
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III. 4. Efficiencies 




4.3 4.4 4.5 
Te„ IMeV] 



FIG. 2. Tcfi spectra for the PMT background events obtained 
from a bifurcated analysis of data from Phase I including the 
best fits to Equation |10| 



rameterized by the following function 
/'pMT(Teff,^i4,/9) =e^^^"■ 



x(e''''+|fe+l|-l) 

•KN{Pli\LUa+UJl{p-QM),Ps 



(10) 



where e, u, b, ujq, uji, and (3^ were parameters determined 
from the fit to the bifurcated data. Af{x\x,(T) represents 
a Gaussian distribution in x with mean x and standard 
deviation a. The uncertainties in the radial parameters 
were obtained from one dimensional scans of the likeli- 
hood function because the magnitude operator distorted 
the likelihood function at 6 = — 1. Compared to the func- 
tion used in the previous analysis [7], wip was replaced 
with ujiip — 0.61) to reduce the correlation between ujq 
and uji . Table |III| lists the PDF parameters from this 
analysis. We observed a weak day/night asymmetry in 
these results, in particular at roughly the la level in the 
Tgff spectrum. Figure [2] shows the fits to the T^tf spec- 
trum for Phase I. 



Table XXI in Appendix [B] lists the uncertainties as- 
sociated with neutron detection. The majority of these 
were identical to previous analyses of these data [71 [H] 
except where indicated below. 

We recently published [M] an analysis based on cal- 
ibration data from Phase III that determined that the 
fraction of neutrons created in NC interactions that were 
detectable with the PMT array, e™'^, was 0.0502 ± 
0.0014. The previous analysis of data from Phase III [8] 
used e™"^ = 0.0485 ± 0.0006, which relied on a Monte 
Carlo based method to determine the uncertainty on the 
neutron detection efficiency. This new analysis avoids the 
dependence on Monte Carlo simulations. 

Similarly, a calibration based method was used to cal- 
culate the fraction of neutrons created by NC interactions 
that were captured in the gas of the NCD array, eJJ'"^^. 
This yielded ej^co ^ o.211±0.005 [H], which had slightly 
better precision than the value of eJJ°° = 0.211 ± 0.007 
used in the previous analysis of data from Phase III [S]. 
We multiplied e^CD ^ correction factor of 0.862±0.004 
in order to determine the efficiency for detecting NC in- 
teractions in the NCD array. The correction factor, av- 
eraged over the duration of Phase III, included the mean 
live fraction of the signal processing systems, threshold 
efficiencies, and signal acceptance due to event selection 
cuts. 

This analysis corrected a 1.2% error in the normaliza- 
tion of the number of NC events observed in the PMT 
array that was in the previous analysis of data from Phase 
III [SI . Because the majority of NC events were observed 
in the NCD array, this normalization error had a rela- 
tively small effect on the measured NC flux. In addition 
we have corrected a 0.1% error in the deuteron density 
used for that analysis. 



TABLE III. PMT background PDF parameters as determined 
by a bifurcated analysis. pi,b is the correlation between the v 
and b parameters. 



Parameter 






Phase 










I-day 


I-night 


Il-day 


II- 


night 


e 


-6.7 


± 1.3 


-5.6 ± 1.0 


-6.3 ± 0.9 


-7.0 


± 0.9 




6.6 


± 0.9 


6.8 ± 1.5 


5.3 ± 1.0 


5.7 


± 1.1 


b 


-1.0 


± 1.3 


3.3 ± 12.0 


-0.3 ± 2.1 


0.5 


± 3.0 


Pub 




0.60 


0.96 


0.91 




0.94 



III. 5. Description of the fit 

The combined fit to all phases was performed using the 
maximum likelihood technique, where the negative log 
likelihood function was minimized using MINUIT [27j . 
The events observed in the PMT and NCD arrays were 
uncorrelated, therefore the negative log likelihood func- 
tion for all data was given by 



The El, spectrum of neutrinos used to derive the 
PDFs for ES, CC, and NC interactions was obtained, in- 
cluding the uncertainties, from Reference [IB]. Our pre- 
vious analysis of data from Phases I and II [7] included 
this uncertainty for the CC and ES PDFs. In this analy- 
sis we also included the effects of this uncertainty on the 
normalization of the NC rate. 



- logLdata = -logipMT - logLNCD, (H) 

where Lpmt and ^ncd, respectively, were the likelihood 
functions for the events observed in the PMT and NCD 
arrays. 

The negative log likelihood function in the PMT array 
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was given by 

- logipMT = Aj (77) 



N 



3 = 1 

npMT 



i=l 



N 



^Aj(rf)/(f,|j,77) 

3 = 1 



,(12) 



where iV was the number of different event classes, ff 
was a vector of "nuisance" parameters associated with 
the systematic uncertainties, Xjijf) was the mean of a 
Poisson distribution for the j"^ class, Xi was the vector 
of event variables for event i, npMT was the total number 
of events in the PMT array during the three phases, and 
f{xi\j, ff) was the PDF for events of type j. 

The PDFs for the signal events were re-weighted based 
on Equations [6] and [7] This was a CPU intensive task 
that was prohibitive for the kernel based PDFs used in 
the previous analysis of data from Phases I and II 7J. 
Therefore, in that analysis, the PDFs were also binned 
based on E^, and these PDFs were then weighted by 
the integral of Equations |6] and [7] within that bin. This 
analysis did not require this approximation when calcu- 
lating the best fit, although an approximation was used 
when "scanning" (described below) the systematic un- 
certainties. As described in Reference [T7], A(7f) was re- 
parameterized such that the Monte Carlo-based PDFs 
did not need to be normalized. 

In the previous analysis of the PMT array data from 
Phase III [^, the nuisance parameters were only propa- 
gated on the PDFs for neutrino interactions, while this 
analysis also propagates these parameters on the PDFs 
for background events. Because the number of back- 
ground events observed in the PMT array was small rel- 
ative to the number of neutrino events, this had a rela- 
tively minor effect on that result. 

The negative log likelihood function in the NCD array 
was given by 



log -Lncd 



O'NCD 



(13) 



where Vjifj) was the mean of a Poisson distribution for 
the j^^ class, uncd was the total number of neutrons 
observed in the NCD array based on the fit described in 
Section p^V.2[ and ctncd was the associated uncertainty. 

The mean number of events for a given class was often 
related to the mean number of events for another class; 
for instance, the number of ES, CC, and NC events in 
each phase were determined from the parameters in Sec- 
tion Umil 

Constraints were placed on various nuisance parame- 
ters and the rate of certain classes of background events. 
Including these constraints, the negative log likelihood 
function was given by 



log L = - log Ldata -t- - (?7 - Ai)'^S" 



\r|-^i), (14) 



where was the value of the nuisance parameters, /2 was 
the constraint on the parameters, and S was the covari- 
ance matrix for the constraints. Appendix [B] lists all of 
the constraints. 

Our previous analysis of data from Phases I and II [7] 
imposed a physical bound on background rates, so that 
they were not allowed to become negative in the fit. 
Without these bounds the background from neutrons 
originating from the AV in Phase II favors a rate whose 
central value was negative, but consistent with zero. The 
uncertainty in the PDFs due to the finite Monte Carlo 
statistics could explain the fitted negative value. The pre- 
viously reported ensemble tests [7] used a central value 
for this background that was more than two statistical 
standard deviations above zero, such that no significant 
effect from the bound was observed. Using a positive 
bound for the backgrounds when ensemble testing with 
Monte Carlo data that does not contain any neutrons 
originating from the AV in Phase II tends to shift 
down on average compared to the flux used to simulate 
the data, as we only obtained background rates that were 
equal to or higher than the values used in the simulations. 
Removing this bound allowed closer agreement between 
the expectation values for the signal parameters and the 
values used in the simulated data. We therefore removed 
these bounds in this analysis to provide a more frequen- 
tist approach that facilitates the combination of the SNO 
results with other experiments. 

In the previous analysis of data from Phases I and II [7] 
the background constraints obtained for the average of 
the day and night rate (e.g. for the ex-situ measurements 
of ^^"'Bi and ^osrpj^-j were applied independently to both 
the day and the night rates, which resulted in a narrower 
constraint on these backgrounds than we intended. This 
analysis correctly applies this as a day and night averaged 
constraint. 

Because of the large number of constraints, and the 
time involved in modifying some PDFs, there were three 
methods for handling the nuisance parameters. Some 
were "floated," i.e. allowed to vary within the MI- 
NUIT minimization of the negative log likelihood func- 
tion. Others were "scanned," where a series of fits were 
performed with different values of the parameter in or- 
der to find the best fit. This process was repeated for 
all scanned nuisance parameters multiple times to con- 
verge on the global minimum of the fit. Finally, some 
were "shifted-and-refitted," where the central values of 
the fit parameters were not varied, but the effect of the 
nuisance parameter on the uncertainties was calculated. 
The method used for each parameter depended on the 
size of that constraints effect on $b and the Ve survival 
probability parameters. Appendix [B] lists how each nui- 
sance parameter was treated. 

In addition to the systematic uncertainties considered 
in previous analyses, this analysis also included a system- 
atic uncertainty due to the finite Monte Carlo statistics 
used to construct the PDFs. We performed 1000 inde- 
pendent fits in which the number of events in each bin of 
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the PDF were drawn from a Poisson distribution. The 
uncertainty due to the finite Monte Carlo statistics was 
determined from the distribution of the fit parameters. 

In order to calculate the total systematic uncertainty 
on the <i>B and the survival probability parameters, we 
applied the shift-and-rcfit 100 times for each parameter 
in order to calculate the asymmetrical likelihood distribu- 
tion for that parameter. We then performed one million 
fits with the nuisance parameters drawn randomly from 
these distributions. The total systematic uncertainty was 
obtained from the resulting distribution of the fit results. 
This is the first time we have applied this procedure, 
which correctly accounts for the combination of asym- 
metrical uncertainties. In order to calculate the effects 
of the day/night or MC systematic uncertainties, respec- 
tively, this procedure was repeated with only the nui- 
sance parameter related to day /night differences or MC 
statistics varied. 

In total the fit included $Bi the five v,^ survival prob- 
ability parameters described in Section III.1[ 36 back- 
ground rate parameters, 35 fioated or scanned nuisance 
parameters, and 82 shift-and-refit nuisance parameters. 

The biases and uncertainties obtained from this anal- 
ysis method were tested using simulated data. The 
number of simulated sets of data was restricted by the 
amount of Monte Carlo data available. For simulated 
data containing neutrino interactions and some back- 
ground classes, we did bias tests with 250 sets of data. 
For simulated data containing neutrino interactions and 
all background classes, we did bias tests with 14 sets of 
data. All of these tests showed the method was unbiased 
and produced uncertainties consistent with frequentist 
statistics. 



III. 6. Crosschecks 

As a crosscheck on the analysis method described 
above, we developed two independent analyses. The first 
crosscheck compared the results from the above method 
run only on data from Phases I and II. This was cross- 
checked against the method described in the previous 
analysis of this data [7] with the most significant changes 
from this analysis included in that previous analysis. The 
results from the two methods were in agreement. 

We developed an alternate Bayesian based analysis 
where the posterior probability distribution was sampled 
using a Markov Chain Monte Carlo (MCMC). This anal- 
ysis was applied to data from Phase III, using the re- 
sults from the maximum likelihood analysis performed 
on data from Phases I and II as a prior. The priors for 
background and neutrino interaction rates had zero prob- 
ability for negative rates and were uniform for positive 
rates. There were two important differences between this 
alternate analysis and the maximum likelihood method 
described above. Firstly, because the systematic uncer- 
tainties were varied in each step of the MCMC, the un- 
certainties included all systematic and statistical uncer- 



tainties. Secondly, this method samples the entire pos- 
terior probability distribution instead of identifying the 
maximum likelihood. Reference [15] provides details of 
this method. As shown in Section |V] the results of the 
Bayesian and maximum likelihood fits agreed. An alter- 
nate Bayesian analysis was also performed with details 
provided in Reference [22] . Though this analysis was not 
used as a direct crosscheck, its results were consistent 
with the analysis presented here. 



IV. NCD ARRAY ANALYSIS 

The NCD array observed neutrons, alphas, and events 
caused by instrumental backgrounds. Because of their 
low stopping power in the gas of the proportional coun- 
ters, electrons rarely triggered the NCD array. A series of 
cuts described in Reference [M] removed the instrumen- 
tal backgrounds. For neutron events i^NCD was peaked 
at approximately 0.76 MeV, with a maximum energy of 
0.85 MeV when including the resolution. i?NCD was less 
than 0.76 MeV if the proton or triton hit the nickel walls 
before losing all their energy. We identified the follow- 
ing two major categories of alpha events: so-called bulk 
alphas, which came from radioactive decays occurring 
throughout the nickel bodies of the proportional counters 
due to the presence of U and Th and their progeny, and 
so-called surface alphas, coming from the decay of ^^"Po, 
which was supported by ^^°Pb that had plated onto the 
inner surface of the nickel bodies. Below 1.4 MeV both 
types of alpha events produced relatively fiat £^ncd spec- 
tra. Due to differences in construction of the strings, the 
number of alpha events varied from string-to-string. 

The previous analysis of data from Phase III [5] distin- 
guished between neutron and alpha events by fitting the 
-2'NCD spectrum. The PDF of i?NCD for neutron events 
was obtained from calibration data, and for alpha events 
it was obtained from simulations. Between 0.4 MeV and 
1.4 MeV the fitted number of alpha and neutron events, 
respectively, were approximately 5560 and 1170. The 
large number of alpha events resulted in both a large sta- 
tistical uncertainty, and a large systematic uncertainty 
due to difficulties in accurately determining the PDF of 
-Encd for alpha events. 

The waveforms of neutron events could be significantly 
broader than those from alpha events, depending on the 
orientation of the proton-triton trajectory. This distinc- 
tion was lessened by the significant tail in the waveforms 
caused by the long ion drift times. In an attempt to 
reduce the number of alpha events, and therefore the un- 
certainties associated with them, we developed four dif- 
ferent particle identification (PID) parameters and a cut 
on these parameters. As described Section p^V.l[ this cut 
reduced the number of events in the strings filled with 
^He (alpha events) by more than 98%, while maintaining 
74.78% of the neutron events. Section |rV.2| describes the 
fit to the -Encd spectrum after this cut. 

These analyses rely heavily on two calibration peri- 
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ods with a Na source distributed uniformly through- 
out the detector [25] , which produced neutrons similar to 
those from *B neutrino NC reactions. These calibrations 
were performed in 2005 and 2006, and were respectively 
referred to as 24Na-2005 and 24Na-2006. A composite 
source of ^"'^ Am and ^Be, referred to as AmBe, produced 
a point-like source of neutrons. This source was posi- 
tioned throughout the detector during six significant cal- 
ibration campaigns spanning Phase III. These data were 
useful for assessing systematic uncertainties associated 
with temporal and spatial variations in the neutron de- 
tection efficiency and PDF of -Encd- 

IV. 1. Particle identification in the NCD array 

Before analyzing the waveforms, the effect of the log- 
arithmic amplifier was removed using parameters deter- 
mined from various calibration pulses in a process re- 
ferred to as de- logging [14] . 

Two particle identification parameters, Pa and pt, were 
based on fitting the waveforms to libraries of known neu- 
tron and alpha waveforms. Each waveform was fitted to 
each library waveform based on a method. The rel- 
ative start time of the event and library waveforms was 
varied to find the minimum x^- In both cases the fits 
did not extend to later times to avoid the effects of ion 
mobility. Both of these particle identification parameters 
were defined by 

log (I). (15) 

where Xa ^^'^ Xrn respectively, were the best x^s from 
the alpha and neutron hypotheses. The libraries used 
to calculate Pa were primarily based on simulation |30j . 
and the was calculated between where the waveform 
first crossed a value equal to 10% of the peak value and 
where it first returned to 30% of the peak value [5T] . 
Figure [3] shows some sample fits. This clearly shows the 
broad waveform for neutrons with a proton-triton trajec- 
tory that was roughly perpendicular to the anode, which 
allows them to be separated from alphas. 

To calculate pb, the neutron library was obtained from 
^^Na-2005 data, and the alpha library was obtained from 
events on the strings filled with '^He [35] • The was 
calculated between where the waveform first crossed a 
value equal to 10% of the peak value and where it first 
returned to 40% of the peak value. The libraries for this 
parameter included events that were used in later studies 
to evaluate performance. We excluded fitting a waveform 
to itself because this would result in a x^ equal to zero, 
i.e. a perfect match. 

The remaining two particle identification parameters, 
Pc and pd, were respectively based on the kurtosis and 
skewness of the waveform after smoothing the waveform 
and deconvolving the effects of ion mobility assuming an 
ion drift time of 16 ns. The skewness and kurtosis were 
calculated using the region between where the waveform 
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FIG. 3. Sample waveforms. The top plot shows a neutron 
waveform (black) obtained from ^*Na calibration data with 
the best fit to the neutron hypothesis (red). The bottom plot 
shows an alpha waveform (black) obtained from a string filled 
with *He with the best fit to the alpha hypothesis (red). The 
vertical lines represent the fit boundaries. 



first crossed a value equal to 20% of the maximum and 
where it first returned to 20% of the peak value. 

Figure [4] shows the distribution of the particle identi- 
fication parameters for known neutron and alpha events. 
The left plot shows that Pa and pi, were highly corre- 
lated, which was unsurprising given their similar defi- 
nitions. This plot also shows that a cut on these two 
parameters (PID cut 1) removes almost all alpha events 
while preserving the majority of neutron events. This 
cut selected events where the alpha hypothesis was sig- 
nificantly worse than the neutron hypothesis. After this 
cut, we recovered approximately 5% of the neutron events 
with a second cut on p^, and pd (PID cut 2) . PID cut 2 
was only applied to events that failed PID cut 1, and 
selected events with high skewness {pd) or low kurtosis 
(pc), he. the waveforms were not symmetric in time or 
had a relatively flat peak. This combined cut, selecting 
events that passed PID cuts 1 or 2, was used for the rest 
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FIG. 4. Distribution of particle identification parameters for 
neutron events (boxes, where the area represents the number 
of events) and alpha events (red marks). The line represents 
the boundary for cuts. PID cut 1 applies to parameters pa 
and pb, and PID cut 2 applies to parameters pc and pd for 
events that failed PID cut 1. 
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FIG. 6. epiD as a function of string for the ^'*Na-2005(black) 
and ^*Na-2006(red) data. Fits to straight lines give 0.7505 ± 
0.0035 with xVNDF of 24.1/29 and 0.7467 ± 0.0018 with 



XVNDF of 49.3/29, respectively. 
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FIG. 5. Bncd spectrum before (brown) and after (red) the 
particle identification cut. From left to right the plots are for 
^^Na calibration data (neutrons), data from strings filled with 
''He (alphas), and data from strings filled with ^He. 



of this analysis. 

Figure [s] shows that the particle identification cut re- 
moves almost all the events on the strings filled with 
"^He, i.e. alpha events, while maintaining the majority 
of the ^"'Na calibration events, i.e. neutron events. This 
also shows that the fraction of alpha events removed by 
the particle identification cut was relatively constant as 
a function of i?NCD- The right most plot of Figure [5] 
shows that the alpha background was significantly re- 
duced, leaving what was clearly mostly neutron events. 

Figure [6] shows the fraction of neutron events surviving 
the combined particle identification cut, epio, as a func- 
tion of neutron capture string for ^''Na-2005 and ^''Na- 



2006 data. Table |IV| shows the average obtained from 
these measurements. The high x^/NDF obtained with 
the ^^Na-2006 data suggests a slight variation in epiD as 
a function of string; however, the correlation between the 
epiD calculated for each string between the 2005 and 2006 
calibrations was only 0.159, which was so small that it 
suggested random string-to-string variation instead of a 
feature of the NCD array. 



TABLE IV. epiD obtained with the ^''Na-2005 and ^''Na-2006 
data. The weighted average included a scaling of the uncer- 
tainty by yp/NDF for the ^''Na-2006 data. 



Data 




X^'/NDF 


^*Na-2005 


0.7505 ± 0.0035 


24.1/29 


^''Na-2006 


0.7467 ± 0.0018 


49.3/29 


Weighted avcrag 


e 0.7478 ± 0.0019 





Table |V] summarizes the systematic uncertainties as- 
sociated with epiD- Based on the methods for deriving 
these systematic uncertainties, we assumed most corre- 
lations were zero. A correlation of 0.50 was assumed 
between the following pairs of systematic uncertainties: 
de-logging and ^^Na uniformity, de-logging and temporal 
variation, pa correction and ^^Na uniformity, Pa correc- 
tion and temporal variation, and Pa correction and de- 
logging. Including these correlations the total absolute 
systematic uncertainty was 0.0065. Combining the sys- 
tematic and statistical uncertainties in quadrature led to 
a total absolute uncertainty of 0.0068. 

The ^^Na calibration data used to calculate epiD had 
a measured variation in the neutron production rate as 
a function of z position of less than 10% between the 
maximum and the value at z = 0. Figure [7] shows that 
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TABLE V. Absolute systematic uncertainties for epm- 



Systematic uncertainty 


'^'^Na uniformity 


0.0010 


Temporal variation 


0.0037 


Contamination 


0.0019 


Do-logging 


0.0018 


Pa correction 


0.0010 


Pf, neutron library 


0.0019 


Total 


0.0065 
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FIG. 7. epiD as function of z for a single string. This was 
typical of all strings. 



the dependence of epiD with source position, as mea- 
sured with the AmBe data, was well approximated by 
a linear function with a maximum deviation compared 
to that at z = of less than 0.01. Combining the possi- 
ble non-uniformity in the ^'^Na source distribution with 
the variation in epiD as a function of z position resulted 
in an absolute systematic uncertainty in epiD due to ^"^Na 
uniformity of 0.0010. The variation due to the x and y 
position non-uniformity was accounted for in the string 
averaging used to calculate epiD- 

The systematic uncertainty in epiD due to temporal 
variations was estimated based on the standard devia- 
tion of epiD calculated from the AmBe data averaged 
over all strings, and calculated at z — assuming a lin- 
ear dependence on z. The systematic uncertainty in epiD 
due to alpha events contaminating the ^^Na calibrations 
was estimated using the number of events with i?NCD be- 
tween 0.9 MeV and 1.4 MeV as an estimate of the alpha 
contamination. The systematic uncertainty in epio due 
to the de-logging process was estimated by recalculating 
epiD with the individual de-logging parameters shifted by 
their estimated uncertainties; because of possible corre- 
lations, the magnitude of the maximum shifts with each 
parameter were added together. 

A correction to the Pa parameter based on ^'^Na and 
AmBe data reduced the spatial and temporal variations 



in this parameter. A systematic uncertainty to account 
for the effect of this correction was estimated by calcu- 
lating epiD assuming a one standard deviation shift in 
the correction to the pa parameter and then combining 
the shifts caused by each string in quadrature, which as- 
sumed that the corrections from string-to-string were not 
correlated. 

The ^^Na-2005 data were used in both the neutron li- 
brary used to calculate pb, and in the determination of 
epiD • Although we did not expect this to bias the calcula- 
tion of cpiD , we conservatively included an additional ab- 
solute uncertainty of 0.0019, half the difference between 
epiD calculated with the 24Na-2005 and 24Na-2006 data. 



IV.2. Method for fitting the NCD array data 

After the particle identification cut, the number of neu- 
tron events was determined from a likelihood fit to a his- 
togram of i?NCD with 50 bins uniformly spaced between 
0.4 and 0.9 McV. 

The PDF of i?NCD for neutron events was obtained 
from ^^Na-2006 data, and for alpha events it was ap- 
proximated by 



Pa{E^cu) =Po 



NCD) 



JVmax 

E 

n=l 



PnPn{E]^CB) 



(16) 



where po and the p„s were fit parameters, P,i(£'NCD) is 
the Legendre polynomial of order n: Pq = 1, Pi — x, 
Pa = 1/2(3x2 - 1), P3 = 1/2(5x3 - 3x), P4 = i/8(35x4 - 
30x2 -h3), .^^ithx = 4(£;NCD[MeV]-0.65). In order to en- 
sure a well defined PDF, negative values of this function 
were set to zero. The fit was repeated with different val- 
ues for the systematic uncertainties associated with the 
Enc d sca le, 

aNCDE^ and resolution, b^'-'^^, (see Equa- 
tion B18 in Appendix |b| selected from a 2-dimension 
scan of these parameters. The point in this 2-dimensional 
scan with the minimum was chosen as the best fit 
point, and the systematic uncertainty associated with 
^NCDE and 5q ^1^^ was obtained from the maximum dif- 
ference in the number of neutron events from the best fit 
point at the Itr contour. 

This fit was performed for values of A^max up to four, 
at which point, based on simulations, the polynomial 
started to fit to fiuctuations in the data. We started 
with the assumption that a zeroth order polynomial was 
a satisfactory fit to the alpha background. If a higher 
order polynomial had a significant improvement in 
then this became the new default, and this was tested 
against higher order polynomials. A significant improve- 
ment in was defined as a decrease in x^ that would 
result in a 32% probability for accepting the higher or- 
der polynomial when the higher order was not a better 
model. This calculation included the fact that testing 
against many different higher order polynomials increases 
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FIG. 8. i^NCD spectrum for events on the strings filled with 
*He after the particle identification cut. The black, red, green, 
and blue lines, respectively, show the PDFs used to simulate 
alpha events for A'^max equal to one, two, three, and four. 



the chances of erroneously choosing a higher order poly- 
nomial, so a larger improvement in was required. Ref- 
erence |33j gives the changes in defined as significant. 
This method was generic to any type of background, 
including instrumental backgrounds, provided they did 
not have features sharper than the assumed background 
shape. 

We tested the bias of this method using simulated data. 
The mean number of neutron events in these sets of simu- 
lated data was based on the number of neutrons obtained 
from the previous analysis of data from Phase III [5] 
and epiD- The i?NCD values for these simulated neutron 
events were obtained from events that passed the particle 
identification cut in the ^**Na-2006 data |j34|- The mean 
number of alpha events in these sets of simulated data 
was based on the number of alphas obtained from the 
previous analysis of data from Phase III [8, and the ap- 
proximate fraction of alpha events removed by the parti- 
cle identification cut. The i^NCD values for the simulated 
alpha events were obtained from events that passed the 
particle identification cut in the strings filled with ^He. 
Because these strings did not have enough events, instead 
of using these events directly, we fitted the limited data 



to polynomials of the form in Equation 16 with A^max 
varied from 1 to 4, and then used these polynomials to 
simulate as many E'ncd values as necessary. In order to 
test extreme possibilities for the alpha event Sncd distri- 
butions, the highest order term from the fit was changed 
by plus and minus 2a, resulting in the eight different 
PDFs shown in Figure [8| The bias was less than 2% for 
all eight alpha PDFs. 

Since only the ^'^Na-2006 data were used to determine 
the PDF of i?NCD for neutrons, we included additional 
systematic uncertainties to account for changes in this 
PDF due to non-uniformity of the ^''Na source and possi- 
ble temporal variations. The size of these systematic un- 



certainties were estimated using AmBe calibration data 
collected at various positions and times to calculate the 
PDF of -Encd for neutron events, and then calculating 
the size of the shift in the reconstructed number of neu- 
tron events. The systematic uncertainties from the tem- 
poral and position variation were summed in quadrature 
to give a total systematic uncertainty of 0.64% on the 
number of neutrons obtained from the fit due to the PDF 
of E'ncd for neutron events. 



RESULTS 



Section |V.1| presents the results from the analysis of 
data from the NCD array in Phase III. Because this was 
a new analysis of this data, we used a statistically-limited 
and randomly-selected one-third subset of the data to de- 
velop the particle identification cut and analysis. Once 
we had finalized all aspects of this analysis we fitted the 
entire set of data from the NCD array in Phase III. After 
completing this full analysis we realized that there was 
an error in the method to calculate the systematic un- 
certainty due to aj^*^^^ and 6^*-^'^^, which was corrected 
in the results presented here. 

The total number of neutron events detected in the 
NCD array obtained from this new analysis of data from 
Phase III was then used as a constraint in the fits to 
the combined data presented in Section |V.2[ The com- 
bined analysis of the three phases also used a statistically- 
limited and randomly-selected one-third subset of the 
data to develop the fitting method. Once we had final- 
ized all aspects of this analysis we fitted the entire set of 
data from all three phases. 



V.l. Results from fit to NCD array data 



Table VI shows the and statistical uncertainty from 
the fit to the i?NCD spectrum for various values of N^ax 



in Equation 16 In general including extra terms in the 
PDF of E'ncd for alpha events should not result in best 
fits with higher but this can occur if the minimiza- 
tion routine finds different local minima. Based on our 
method for choosing the value of A^max representing the 
point where improvements in fit quality cease, the best fit 
occurs when A'^ax — 4. This corresponded to the maxi- 
mum value of A^max considered before performing the fit, 
so to check that larger values of A'^max did not produce 
better fits, we also fitted with A'max equal to five and six, 
as shown in Table |VI| These fits did not produce better 
results. 

Figure [9] shows the best fit of the i?NCD spectrum. Al- 
though the best fit turns down at higher values of -Encd 
the parameters were consistent with a flat PDF in that 
region. This variation in the allowed PDF was reflected 
in the increased statistical uncertainty with large A^max- 

For the fit with A^max = 4 the systematic uncertainty 
due to ^'^^ and 6^^^^ was 5 neutrons. Combining this 
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TABLE VI. Bind la statistical uncertainty for various val- 
ues of A'^max in Equation 1 16[ 



A^max X^/NDF Stat, uncertainty 



54.92/48 

1 56.72/47 
47.63/46 
41.78/45 
40.20/44 
40.34/43 
40.41/42 



4.2% 
4.2% 
5.5% 
6.5% 
6.9% 
9.4% 
9.2% 




3.4 0.45 0.5 0.55 0.6 0.65 0.7 0.75 0.8 0.85 

Ener£ 



0.9 
,y [MeV] 



FIG. 9. The fitted Bncd spectrum after the particle identi- 
fication cut. The thick black line is the best fit. The blue 
and red lines are the best fitted neutron and alpha spectra, 
respectively. 



with the systematic uncertainty in the PDF of £^ncd for 
neutrons, the statistical uncertainty in the fit, and divid- 
ing by epiD , the total number of neutrons observed in the 
NCD array equals 1115 ± 79. The previous analysis of 
data from Phase III gave 1168 neutrons in the NCD array, 
with similar uncertainty [S]. That analysis had a large 
background due to alpha events, which made the assess- 
ment of the systematic uncertainty on the fitted number 
of events more challenging. The result presented here 
avoids that problem by eliminating most of the back- 
ground from alpha events and allowing a very general 
PDF to describe the i?NCD spectrum for any remaining 
background events. Since the particle identification cut 
removed almost all alpha events, the fitted number of 
neutron events had a small to moderate correlation with 
the previous analysis of this data. 



V.2. Results from combined fit to all data 

For the combined fit to all data using the maxi- 
mum likelihood technique. Table [VlT| shows $b and the 
Ve survival probability parameters as defined in Equa- 
tions m and [7] of Section HILTI Table IVIIll shows the 



correlation between these parameters. The combined 
fit to all data from SNO yielded a total flux of ac- 
tive neutrino flavors from ®B decays in the Sun of 
$B=(5.25±0.16(stat.)+°;J3(syst.)) x lO^cm-^s-i. Dur- 
ing the day the v^, survival probability at 10 MeV was 
Co = 0.317 ± 0.016(stat.) ± 0.009{syst.), which was in- 
consistent with the null hypothesis that there were no 
neutrino oscillations at very high significance. Using the 
covariance matrix obtained from this combined analysis 
we can compare the best fit to various null hypotheses. 
The null hypothesis that there were no spectral distor- 
tions of the Ve survival probability (i.e. ci — 0, C2 — 0, 
ao = 0, ai = 0), yielded Ax^ = 1.97 (26% C.L) compared 
to the best fit. The null hypothesis that there were no 
day/night distortions of the survival probability (i.e. 
ao = 0, ai = 0), yielded Ax^ = 1.87 (61% C.L.) com- 
pared to the best fit. 



TABLE VII. Results from the maximum likelihood fit. Note 
that "^B is in units of x 10® cm~^s~^. The D/N systematic un- 
certainties includes the effect of all nuisance parameters that 
were applied differently between day and night. The MC sys- 
tematic uncertainties includes the effect of varying the number 
of events in the Monte Carlo based on Poisson statistics. The 
basic systematic uncertainties include the effects of all other 
nuisance parameters. 



Best fit 


Stat. 


Systematic uncertainty 






Basic D/N MC Total 



Co 

ci 

C2 

ao 
ai 



5.25 ±0.16 
0.317 ±0.016 
0039 +0 0065 

U.UUda _0.0067 

-0.0010±0.0029 
0.046 ±0.031 
-0.016 ±0.025 



+ 0.11 
-0.12 

+0.008 

-0.010 
+ 0.0047 +0.0012 
-0.0038 -0.0018 
+0.0013 +0.0002 
-0.0016 

+0.007 

-0.005 

+0.003 

-0.006 



±0.01 
±0.002 



0.0003 

±0.012 



+ 0.01 
-0.03 
+ 0.002 
-0.001 
+ 0.0004 
-0.0008 
+ 0.0004 
-0.0002 
+ 0.002 
-0.003 



±0.009 ±0.002 



+0.11 
-0.13 

±0.009 
±0.0045 

+ 0.0014 
-0.0016 
+0.014 
-0.013 
+ 0.010 
-0.011 



TABLE VIII. Correlation matrix from the maximum likeli- 
hood fit. 

*B Co ^^1 ^^2 Oq Ql 

$B 1.000 -0.723 0.302 -0.168 0.028 -0.012 

Co -0.723 1.000 -0.299 -0.366 -0.376 0.129 

ci 0.302 -0.299 1.000 -0.206 0.219 -0.677 

C2 -0.168 -0.366 -0.206 1.000 0.008 -0.035 

ao 0.028 -0.376 0.219 0.008 1.000 -0.297 

ai -0.012 0.129 -0.677 -0.035 -0.297 1.000 



Figure 10 shows the RMS spread in P^^{E,^) and 
Aee{E^), taking into account the parameter uncertain- 
ties and correlations. This also shows that the maximum 
likelihood analysis was consistent with the alternative 
Bayesian analysis. Reference [3S] contains all steps of 
the MCMC fit after the fit had converged. 

Figures [Tl] [T2j and [13) respectively, show one- 
dimensional projections of the fit for Phase I, II, and 
III. 
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r^igh,/«bins: 17.109/21 x;?,/«bi„s: 19.669/21 



FIG. 10. RMS spread in Pe,i{E^) and AeeiE^), taking into 
account the parameter uncertainties and correlations. The 
red band represents the results from the maximum likelihood 
fit, and the blue band represents the results from the Bayesian 
fit. The red and blue solid lines, respectively, are the best fits 
from the maximum likelihood and Bayesian fits. 



VI. NEUTRINO OSCILLATIONS 

Solar models predict the fluxes of neutrinos at the 
Earth, but the flavors of those neutrinos depends on the 
neutrino oscillation parameters and the details of where 
the neutrinos were produced in the Sun. Section [VLT] de- 
scribes how the flavor components of the neutrino fluxes 
were calculated, and Section |VI.2| describes how these 
predictions were compared to results for ^q, 
and Af.e(E^) presented here, and with other solar neu- 
trino experiments. Reference [36] provides further details 
on the neutrino oscillation analysis presented here. 

We considered two different neutrino oscillation hy- 
potheses in this analysis. For comparison with previous 
SNO analyses. Section |VI.3| presents the so-called two- 
flavor neutrino oscillations, which assumed ^13 = and 
had two free neut rino oscillation parameters, 9i2, and 
Am2i. In Section VI. 4 we also considered the so-called 
three-flavor neutrino oscillations, which had the follow- 



ing three free neutrino oscillation parameters: 612, ^13, 
and Amjj. Note that the mixing angle, 623, and the 
CP-violating phase, S, are irrelevant for the neutrino os- 
cillation analysis of solar neutrino data. The solar neu- 
trino data considered here was insensitive to the exact 
value Am|]^, so we used a fixed value of ±2.45 x 10'^ eV^ 
obtained from long-baseline accelerator experiments and 
atmospheric neutrino experiments |37| . 
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FIG. 11. Projection of the T^s, p, cosOq, and /3i4 for the 
Phase I data. Day events hollow circles and dashed lines. 
Night events filled circle and solid lines. Note that the sharp 
break in the data in the top panel at 5 MeV arises from change 
of bin width. 



VI. 1. Solar neutrino predictions 

Predicting the solar neutrino flux and spectrum 
for all neutrino flavors requires a model of the neutrino 
production rates as a function of location within the 
Sun, and a model of the neutrino survival probabili- 
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FIG. 12. Projection of the Tofi, p, cosOq, and /3i4 for the 
Phase II data. Day events hollow circles and dashed lines. 
Night events filled circle and solid lines. Note that the sharp 
break in the data in the top panel at 5 MeV arises from change 
of bin width. 



ties as they propagate through the Sun, travel to the 
Earth, and then propagate through the Earth. For con- 
sistency with previous calculations, and because of the 
conservative model uncertainties, we used the BS05(OP) 
model |21j to predict the solar neutrino production rate 
within the Sun rather than the more recent BPS09(GS) 
or BPS09(AGSS09) models 6. Reference 05j provides 
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FIG. 13. Projection of the T^a, p, and cos 6© for the data 
from Phase III. Day events hollow circles and dashed lines. 
Night events filled circle and solid lines. 



the results presented below assuming these other solar 
models. We used the E^, spectrum for neutrinos from 
Reference [26j, and all other spectra were from Refer- 
ence [35]. 

In previous analyses we used numerical calculations to 
construct a lookup table of neutrino survival probabihty 
as a function of the neutrino oscillation parameters. Such 
a table was still used to study the entire region of neu- 
trino oscillation parameters. Previous analyses of SNO 
data combined with other solar neutrino experiments left 
only the region referred to as LMA. This analysis used an 
adiabatic approximation when calculating neutrino oscil- 
lation parameters in that region. We verified that these 
two calculations gave equivalent results for a fixed set of 
neutrino oscillation parameters in the LMA region. Due 
to the improved speed of the adiabatic calculation we 
could scan discrete values of both Am^i and E^, whereas 
the lookup table used previously was calculated at dis- 
crete values of Am^i/ E^, which resulted in small but 
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observable discontinuities. 

We also updated the electron density as a function of 
Earth radius, which affects the survival probability of 
neutrinos propagating through the Earth. Previous SNO 
analyses used the PREM [39, model, which averages over 
the continental and oceanic crust. When neutrinos en- 
ter the SNO detector from above they must have prop- 
agated through continental crust; therefore, we modified 
the Earth density profile to use PEM-C [40], which as- 
sumes continental crust for the outer most layer of the 
Earth. Because this significantly affected neutrinos only 
during a short period of each day, this had a negligible 
effect on the calculated neutrino survival probability. 



VI. 2. Analysis of solar neutrino and KamLAND 
data 



To compare the survival probability parameters cal- 
culated in Section |V] with a neutrino oscillation predic- 
tion it was necessary to account for the sensitivity of the 
SNO detector. Table |XXIl| in Appendix [C] gives S{E^), 
the predicted spectrum of detectable by the SNO de- 
tector after including all effects such as the energy de- 
pendence of the cross-sections, reaction thresholds, and 
analysis cuts, but not including neutrino oscillations. We 
multiplied S{E^) by the predicted neutrino oscillation 
hypothesis distortions, and fitted the resulting spectrum 
to S{E^) distorted by Equations [6] and [t] We then calcu- 
lated the between the results from this fit and our fit 
to the SNO data presented in Section [V.2| This calcula- 
tion used the uncertainties and correlation matrix from 
the fit to SNO data, but did not include the uncertain- 
ties from the fit to the distorted S{E^) as this does not 
represent a measurement uncertainty. 

The ■^^s calculated as a function of the neutrino 
oscillation parameters. The best fit was determined from 
the parameters resulting in the minimum x^, and the 
uncertainties were calculated from the change in from 
this minimum. Tests with simulated data revealed that 
this method produced neutrino oscillation parameters 
that were unbiased and produced uncertainties consis- 
tent with frequentist statistics. 

The following additional solar neutrino results were 
used in calculating the results from solar neutrino ex- 
periments: the solar neutrino rates measured in Ga |41j . 
and CI [H] experiments, the rate of ^Be solar neu- 
trinos measured in Borexino [43], and the rates and 
recoil electron energy spectra of ^B neutrino ES re- 
actions measured in Borexino ^44j, SuperKamiokande- 
I [IS], SuperKamiokande-II [H], and SuperKamiokande- 
III [47] . The last two SuperKamiokande results were split 
into day and night, and the first SuperKamiokande result 
was split into multiple periods of the day. The difference 
in the day and night rate of ''Be solar neutrinos recently 
measured in Borexino |48j . and the recent measurement 
of the *B neutrino spectrum in KamLAND were not 
included, but these results would not significantly change 
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FIG. 14. Two-flavor neutrino osciUation analysis contour us- 
ing only SNO data. 



the results reported here. For a given set of neutrino os- 
cillation parameters and 4>b, the predictions for the set 
of experiments were calculated and compared to their 
results. This comparison was added to the described 
above, and the resulting was then minimized with re- 
spect to $6- The same procedure as above was used to 
determine the best fit values and uncertainties. 

The KamLAND experiment observed neutrino oscil- 
lations in DgS from nuclear reactors. By assuming 
CPT invariance we can directly compare these results 
with the neutrino oscillations observed with solar neu- 
trinos. Because this was a completely independent re- 
sult, the lookup table of x^ as a function of 0i2, 0i3, and 
Am^i published by the KamLAND collaboration [SU] was 
added directly to the values calculated from the solar 
neutrino analysis, and the same procedure was used to 
determine the best fit values and uncertainties. 



VI. 3. Two-flavor neutrino oscillation analysis 

Figure [M] shows the allowed regions of the 
(tan^ 6i2, Am2i) parameter space obtained with the re- 
sults in Tables IVIII and IVIIII SNO data alone cannot 
distinguish between the upper (LMA) region, and the 
lower (LOW) region, although it slightly favors the LMA 
region. 

Figure [TSj shows the allowed regions of the 
(tan^ 6*12, Ato2i) parameter space obtained when the 
SNO results were combined with the other solar neu- 
trino experimental results, and when this combined solar 
neutrino result was combined with the results from the 
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FIG. 15. Two-flavor neutrino osciUation analysis contour us- 
ing both solar neutrino and KamLAND (KL) results. 



KamLAND (KL) experiment. The combination of the 
SNO results with the other solar neutrino experimental 
results eliminates the LOW region, and eliminates the 
higher values of Am^^ in the LMA region. 

Table IIXI summarizes the results from these two- flavor 
neutrino analyses. 



TABLE IX. Best-fit neutrino oscillation parameters from a 
two-flavor neutrino osciUation analysis. Uncertainties listed 
are ±lcr after the was minimized with respect to all other 
parameters. 
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VI. 4. Three-flavor neutrino oscillation analysis 

Figure [16] shows the allowed regions of the 
(tan^ 012, Amji) and (tan^ 6*12, sin^ ^13) parameter spaces 
obtained from the results of all solar neutrino exper- 
iments. It also shows the result of these experiments 
combined with the results of the KamLAND experiment. 
Compared to the result in Figure [15] this clearly shows 
that allowing non-zero values of ^13 brings the solar neu- 
trino experimental results into better agreement with the 
results from the KamLAND experiment. 

Figure [17] shows the projection of these results onto 
the individual oscillation parameters. This result shows 



FIG. 16. Three-flavor neutrino oscillation analysis contour 
using both solar neutrino and KamLAND (KL) results. 



that due to the different dependence between tan^ Q\i 
and sin^ ^13 for the solar neutrino experimental results 
and the KamLAND experimental results, the combined 
constraint on sin^ ^13 was significantly better than the 
individual constraints. 

Table 1x1 summarizes the results from these three- flavor 
neutrino oscillation analyses. Tests with the inverted hi- 
erarchy, i.e. negative values of Atoji, gave essentially 
identical results 
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TABLE X. Best-fit neutrino oscillation parameters from 
a three-flavor neutrino oscillation analysis. Uncertainties 
listed are ±lcr after the ^^is minimized with respect 
to all other parameters. The global analysis includes So- 
lar-HKL-hATM-HLBL-hCHOOZ. 



tan 



Analysis 
Solar 

Solar-I-KL 0.446 
Global 



l[cV^] 

0.436t-,-« - X 10- 

^Voll 7.4i;«;- X 10- 



An 



sin^ ei3(xl0-^) 
< 5.8 (95% C.L.) 



2 5 + 1-** 



< 5.3 (95% C.L.) 

,+0.88 
--0.55 



2. 021 



1.00 
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a significance of approximately 2.7a. A combined anal- 
ysis of all LBL and atmospheric (ATM) results, and the 
results from the CHOOZ j53| experiment was performed 
by Fogli et al. [M]. Because the LBL-hATM+CHOOZ 
analysis was insensitive to 612, and because the solar 
neutrino-|-KamLAND analysis was insensitive to Am^^ 
we can simply add their projections of Ax^ onto O13. 
Table |X] and Figure [19] show the results of that combi- 
nation. This shows that the LBL-hATM-|-CHOOZ ex- 
periments currently have better sensitivity to 613 than 
the combined solar and KamLAND experiments, but the 
combination of all experiments gives a slightly improved 
determination of ^13, hinting at a non-zero value. 



VII. DISCUSSION AND CONCLUSIONS 




FIG. 17. Projections of the three-flavor neutrino oscillation 
parameters determined from Figure [16] The horizontal lines 
representing the Ax^ for a particular confidence level, are for 
two-sided intervals in plot (a) and (b), and one-sided intervals 
in plot (c). 



Figure 18 shows the measured solar i/e survival prob- 
ability as a function of Ei,. At higher E^, the results of 
this analysis provide the best constraints on the survival 
probability. All solar results are consistent with the LMA 
neutrino oscillation hypothesis. 

Recent results from the T2K M] and MINOS (55] long- 
baseline (LBL) experiments indicate a non-zero ^13 with 



By developing a particle identification cut to analyze 
data from Phase III of the SNO experiment we measured 
1115 ± 79 neutrons. Eliminating most of the background 
from alpha events and allowing a very general PDF to de- 
scribe the i?NCD spectrum for any remaining background 
events made this analysis less sensitive to background un- 
certainties than our previous analysis of these data. 

Combining data from all phases of the SNO experiment 
we measured a total flux of active flavor neutrinos from 
^B decays in the Sun of (5.25 ± 0.16(stat.)^Q :[3(syst.)) x 
10^ cm~^s~^. We improved the handling of a number of 
systematic uncertainties in this analysis compared with 
our previous analyses of these data. This result was con- 
sistent with but more precise than both the BPS09(GS), 
(5.88±0.65) X 10^ cm-2s-\ and BPS09(AGSS09), (4.85± 
0.58) X lO^cm^^s^^, solar model predictions [HI- 

The precision of the Ve survival probability parame- 
ters was improved by approximately 20% compared to 
our previously reported results due to the additional 
constraint provided by the data from Phase III. Dur- 
ing the day the i^e survival probability at lOMeV was 
Co = 0.317 ± 0.016(stat.) ± 0.009(syst.), which was in- 
consistent with the null hypothesis that there were no 
neutrino oscillations at very high significance. The null 
hypotheses that there were no spectral distortions of the 
i^e survival probability (i.e. ci = 0, C2 = 0, aq = 0, 
oi — 0), and that there were no day/night distortions of 
the survival probability (i.e. ag — 0, ai = 0) could not 
be rejected at the 95% C.L. 
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FIG. 18. Various solar survival probability measurements 
compared to the LMA prediction for *B neutrino. Using the 
results from Section [yi] of this paper, the dashed line is the 
best fit LMA solution for ®B neutrinos and the gray shaded 
band is the la uncertainty. The corresponding bands for i^eS 
from the pp and ^Be reactions (not shown) are almost iden- 
tical in the region of those measurements. The blue shaded 
band is the result of the measurement the ®B neutrino sur- 
vival probability reported here. The red point is the result of 
the Borexino measurement [43j of the survival probability for 
Ve.s produced by '^Be-l-e" — >■ ^Li+Ve reactions in the Sun. The 
blue point is the result of various measurements [41] of the 
survival probability for v^s produced by p -|-p — >■ -|- e+ -I- i^e 
reactions in the Sun; note that these measurements did not 
exclusively measure this reaction, so the contribution from 
other reactions were removed assuming the best fit LMA so- 
lution, and so actually depends on all solar neutrino results. 
The uncertainty in absolute flux of the subtracted reactions 
was included in the calculation of the total uncertainty of 
this point, but the uncertainty due to the neutrino oscillation 
probability of these reactions was not. The uncertainty due to 
the normalization of the two points by the expected flux was 
included. For clarity, this plot illustrates the LMA solution 
relative to only a subset of the solar neutrino experimental 
results. 



two-flavor neutrino oscillation analysis yielded 
= (5.6^^;!^) X 10-5 and tan^ 6*12 = 0A27to°w- 
A three-flavor neutrino oscillation analysis combining 
this result with results of all other solar neutrino exper- 



iments and the KamLAND experiment yielded Arn^i = 
(7.4ltEJ;2i) X 10-5 eV2, tan^ 6*12 = 0.446t°;°^°, and 
^ ^ = {2.5tli) X 10"^- This implied an upper bound 
ha < 0.053 at the 95% C.L. 
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TABLE XIL *B neutrino interactions scaling factors. 



Interaction Day/Night Scaling factor 



CC.ESe 

CC.ESc 

ES^T- 

NC 



Day 
Day 
Night 
Niglit 



Day-I-Niglit 



*b'[1 - PL{E.)' 
*B'P°e(-B^)' 
*B'fl - P^iEu)' 



(1 — Pes{Ev)) and our measurement of Pf^{Ejj) vi^ould be 
scaled by 1/(1 -Pes 

We note that the approximations made in this analysis 
were also valid for effects involving sterile neutrinos with 
spectral distortion that are significant only below about 
4 MeV and with small day-night effects for the NC detec- 
tion process. In this case the principal additional effect 
would be a further renormalization of the NC interaction 
rate. 



Appendix B: Constraints on the likelihood fit 

Systematic uncertainties were propagated on the PDFs 
as described in this section. 



Appendix A: Sterile neutrinos 

If we assume a sterile neutrino, where the probability 
of an electron neutrino oscillating into a sterile neutrino, 
PesiEy), was the same during the day and night, then 
the scaling factors given in Table |l] of Section |III.1| are 
replaced with those in Table |XI| 



TABLE XL *B neutrino interactions scaling factors including 
a probability of an electron neutrino oscillating into a ster- 
ile neutrino, which was the same during the day and night. 



Teff PDFs 



Table XIII lists the systematic uncertainties on the 
PDFs of Teff. For Phases I and II the modified Tog, Tog', 
during the day was obtained from 



Teff' = (1 + ^ 



'0 c 



I^-Af„/2-Af„/2)reff. (Bl) 



2 + Aee(fi„) 



The signs of the ylf„ and Af^^ terms were reversed for 
night. The uncertainty due to Tgfj resolution was ob- 
tained by convolving Teff' with Gaussians centered at 



spectrum of -B^ detectable by the SNO detector after includ 
ing the energy dependence of the cross-section 

Interaction Day/Night Scaling factor 



and f{E,) was the predicted ^ero with widths of , and crf„ and crf„ scaled by 



CC.ESe 

ES^T- 
CC, ESc 

ES^-r 

NC 



Day 
Day 
Night 
Night 

Day+Night *b 



*B-P,,(B,.) 

*b[1 - Pee(-E„) - Pcs{E^)] 

*BP"e(-E,.) 
*B[1-P"e(-E.)- Pea 



If Pes{E^) was a constant as a function of E^, and 
defining $b' = $6(1 - Pes{E,)), and PHE,)' = 
" we obtain the scaling factors given in Ta- 



l~ P^s(E ) 

bielxni 



Notice that scaling factors in Table |XII| are equiv- 
alent to those in Table |T] except our measurement of 
the neutrino flux would be the true flux scaled by 



parameterized detector resolution. This resolution was 
applied first to just day events and then to just night 
events. Differences in the neutrino parameters between 
the shifted fits and the central fit were taken as the re- 
sulting uncertainties. 

During Phase III the modified Teg, Teff", was obtained 
from 

Teff" = Teff' + 5f (1 - S|j2)(Teff' - Tg), (B2) 

where Tg was the true MC energy. For the NC and ex- 
ternal neutrons PDFs Tg was constant and equal to the 
mean fitter energy 5.65 MeV. Teff', was obtained from 

Teff' = (1 -t- af , + )(1 - Af„/2 - A|,/2)Teff ,(B3) 

where the signs of the A^^ and A^^^ terms were reversed 
for night. 
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TABLE XIII. Tcti PDF systematic uncertainties. 



Parameter Description 



Events Phase Nominal Variation 



Fit Value Application 



0004+" ""'"' 
-0 0007+" ""^** 

U.UUU ( _Q QQ3Q 

0001+" ""^'' 

U.UUU1_Q QQ27 
0.0065_(, oog4 

N/A 
N/A 
N/A 
N/A 
N/A 
N/A 
N/A 

n ooo»;+"-""^^ 

n nn'?«+" ""ss 
2l4+" "23 

"■^"■^-0.041 
+ 0.017 
U.100_o 019 

0109+"-"^"^ 

U.UlU»_o.0100 

N/A 
N/A 
N/A 
N/A 
N/A 
N/A 
N/A 
0.000 ± 0.012 



A', 
Al 



[MeV] 
[MeV] 
[MeV] 



Tcff scale 
Tcff scale 
Tcff scale 
Teff scale 

Teff scale non-linearity with Tcff 
Teff scale diurnal variation 
Toff scale directional variation 
Teff scale directional variation 
Teff scale diurnal variation 
Tcff scale directional variation 
Teff scale directional variation 
Tcff scale diurnal variation 
Tcff scale directional variation 
Tcff resolution 
Tcff resolution 
Teff resolution 
Teff resolution 
Teff resolution 

Teff resolution diurnal variation 
Tcff resolution directional variation 
Tcff resolution directional variation 
Tcff resolution diurnal variation 
Tcff resolution directional variation 
Teff resolution directional variation 
Tcff resolution directional variation 



all 
all 
all 
all 
all 
all 
CC 
ES 
all 
CC 
ES 
all 
ES 
all 



all 
CC 
ES 
all 
CC 
ES 
ES 



I, II, 
I 

II 

III 

I, II, 

I 

I 

I 

II 

II 

II 

III 

III 

I 

II 

II 

III 

III 

I 

I 

I 

II 
II 
II 
III 



III 



III 

















0.155 
0.168 
0.154 
0.0119 
0.016184 










±0.0041 

+ 0.0039 
-0.0047 
+ 0.0034 
-0.0032 
±0.0081 

±0.0069 
±0.0032 

±o.ooof 

±0.009;^ 
±0.004 
±0.000^ 
±0.007!^ 
±0.0038 
±0.0099 

+0.041 
-O.OSO 
+ 0.041 
-O.OSO 

±0.018 
±0.010 
±0.014 

±0.00'3 
±0.001^ 
±0.006'* 

±0.005 

±o.ooi:F 

±0.0K^ 
±0.012 



shift- 
shift- 
shift- 

shift- 
shift- 



shift 
shift 



shift 
shift 



scanned 
scanned 
scanned 
scanned 
and- refit 
and- refit 
and-rcfit 
N/A 
and- refit 
and-rcfit 
N/A 
scanned 
scanned 
scanned 
scanned 
scanned 
scanned 
N/A 
and-rcfit 
and-rcfit 
N/A 
and-rcfit 
and-rcfit 
N/A 



^ Correlation of -1. 
^ Correlation of -1. 

Correlation of 1. 

Correlation of -1. 

Correlation of -1. 



2. /3i4 PDFs 



3. p PDFs 



Table XIV lists the systematic uncertainties on the 
PDFs of /3i4. During the day these uncertainties were 
applied using 



PiA = /?i4(l + a^^^ + eg (Teff - 5.589[MeV])) 



(B4) 



The signs of the A^^^ and v4^-^ terms were reversed for 
night. For electrons and 7- rays the systematic uncertain- 
ties associated with resolution were applied as 



(B5) 



Table XV lists the systematic uncertainties on the 
PDFs of p. The radius was modified by 



[{x"[ 



+ (2/"[cni])2 + (z"[cni])' 



i3/2 



6003 



(B6) 



where x", y", and z" were the modified cartesian coordi- 
nates, as described below. Each event was weighted by a 
factor l + cP X (Teff - 5.05[MeV]). 

During Phases I and II, x" , y", and z", respectively, 
were obtained by convolving x', y', and z' with Gaussians 
centered at zero with widths of , cr^, cr^, (t^„, or tr^-^. 
These resolutions were applied first to only day events 
and then to only night events, a;', y', and z' were obtained 
from 



where Pu" was the value of including all of the sys- 
tematic uncertainties, and /3i4 a parameterized average 
value of /3i4 for the PDF. For neutrons in Phase II the 
systematic uncertainties associated with resolution were 
applied as a convolution with a Gaussian centered at zero 
with a width of a^^'^ . This correction can only be applied 
in the positive direction, and the negative fit uncertain- 
ties were inferred to be the same as the positive ones. 



x'^il + a1~ AP,J2 - A2„/2) (a; + (B7) 

y' = {l + a1- A^Jl - A'',j2){y + al) (B8) 
z'^{l + a{ + al-A^,j2-AP,j2) 

(z + ag). (B9) 

The signs of the ^^^^ and A^^-^ terms were reversed for 
night. 
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TABLE XIV. Pi4 PDF systematic uncertainties. 



Paramete 
~P 



4 

4 
4 



Description 



Events Phase Nominal Variation 



Fit Value Applieation 



/3l4 


scale 


all 


I 


/3l4 


scale 


e^ 


II 


/3l4 


scale 


n 


II 


/3l4 


scale non-linearity with Tcff 


all 


I, II 


/3l4 


offset diurnal variation 


all 


I 


014 


offset directional variation 


CC 


I 


01i 


offset directional variation 


ES 


I 


/3l4 


offset diurnal variation 


all 


II 


/3l4 


offset directional variation 


CC 


II 


/3l4 


offset directional variation 


ES 


II 


/3l4 


resolution 


all 


I 


/^14 


resolution 


e^ 


II 


;9i4 


resolution 


n 


II 



-0.0081 


-0.0144 
0.00275597 









0.0150 



±0.004: 

±0.002^ 



N/A 



00102+° '"'"^ 

+ 0.00.^^ _ 
-0.0022 ^■^-^'-"^-0.0025 

±0.00069 0.0020lto:51™44 

±0.0043 
±0.0003t^ 
±0.003'*^ 

±0.0043 
±0.0003^ 
±0.003'^ 
±0.004;^ 
±0.005''^ 

±0.0045 



N/A 
N/A 
N/A 
N/A 
N/A 
N/A 
N/A 
N/A 
N/A 



scanned 
N/A 
scanned 
scanned 
shift-and-refit 
sliift-and-refit 
N/A 

shift-and-refit 
shift-and-refit 
N/A 

shift-and-refit 
N/A 

shift-and-refit 



^ Correlation of 1. 
^ Correlation of -1. 

Correlation of -1. 

Correlation of 1. 



During Phase III the uncertainties were apphed as 

x" = x' + {b^^y + bl'^z + hl'yz^){x' - xg) (BIO) 
y" = y' + + br-^ + brz^){y' - y,) (Bll) 
z" ^ z' + {bl+blz)(z' - zg), (B12) 

where Xg, yg, Zg were the true MC positions, x', y', and 
z' were obtained from 

x' = {l + - - A,u,/2){x + ag) (B13) 

y' = (1 + 0(1 - - A,u,/2){y + al) (B14) 

z' = (1 + < + aj)(l - Adn/2 - Ad,r/2) 

(z + flo)- (B15) 

The signs of the A'^^ and A^^^ terms were reversed for 
night. 



tained from 

£^ncd' = (l + af~°°)i^NCD. (B18) 

For each event in the PDF one hundred random numbers 
drawn from a Gaussian centered at zero with a width of 
bg^'^" E^cu were used to construct a new PDF. 

6. Background constraints 

Table pCVIIII gives the constraints on the backgrounds 
in Phases I and II. Table |XIX| gives the constraints on 
the backgrounds in Phase III. 

7. PMT background PDF 



cos 00 PDFs 



Table |XVI| lists the systematic uncertainties on the 
PDFs of cos 00 . For Phases I and II the modified cos 9q , 
COS00', was obtained from 

cos 6*0' = 1 + (1 + - yl^,,j2)(cos6'o - 1). (B16) 

For Phase III the modified cos 610, cos 6*0', was obtained 
from 

COS00' = 1 + (1 + a^)(l - Alj2){cos0Q - 1XB17) 

If the transformation moved cos 610' outside the range 
[— 1, 1], COS00' was given a random value within this in- 
terval. 



PDFs 



Table |XVII| lists the systematic uncertainties on the 
PDFs of i?NCD- The modified E'ncd, -Encd', was ob- 



Table PCX| shows the constraints on the analytical PDF 
given by Equation [lO] in Section [lIL3l for PMT back- 
ground events. 



8. Neutron detection efficiencies 



Table IXXll shows the constraints on the neutron detec- 
tion efficiencies. NC interactions observed with the PMT 
array in Phases I and II were weighted by 



1 -L pPMT I PMT 
' n ' n.corr. 



(B19) 



NC interactions observed with the PMT array in Phase 
III were weighted by 



1 + e; 



PMT 



(B20) 



NC interactions observed with the NCD array in Phase 
III were weighted by 



1 + e 



NCD 



(B21) 
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TABLE XV. p PDF systematic uncertainties. 



Parameter Description 



Events Phase 



Nominal 



Variation 



+ ().U01() 

-0.0057 

+ 0.0040 

-0.0 

+ 0.0004 

-0.0034 

+ 0.0003 

-0.0025 

+ 0.0029 

-0.0077 

+ 0.0015 

-0.0012 

+ 0.0085 

-0.0049 

+ 0.0041 

-0.0048 

+ 0.0088 

-0.0067 

±0.002 

±0.000.= 

±o.oo; 

±0.0(fe' 
±0.000 
±0.001 
±0.00 Is* 
±0.0018 

+ 1.15 
-0.13 
+2.87 
-0.17 
+2.58 
-0.15 
+0.62 
-0.07 
+2.29 
-0.09 
+3.11 
-0.16 
±4.0 

±4.0 
±4.0 
±3.3 
±2.2 
±1.5 
±3.1 
±3.4 
±5.3 
±0.02; 
±6.1 X 10 
±2.0 X 10 
±0.02J 
±0.83 X 10 

±6.37* 
±7.21 

±i.of 

±3.3e' 



Fit Value Application 



4. 

At 



[cm 
[cm 
[cm 
[cm 
[cm 
[cm 
[cm 
[cm 
[cm 
[cm 
[cm 
[cm 
[cm 
[cm 
[cm 



" [cm 
^ [cm~ 

[cm-i] 
n [cm] 

ir [C"^] 

[cml 



p scale all I 

z scale all I 

p scale all II 

z scale all II 

p scale all III 

z scale all III 

p scale non-linerity with Toff all I 

p scale non-linerity with T^ff all II 

p scale non-linerity with T^ff all III 

p scale diurnal variation all I 

p scale directional variation CC I 

p scale directional variation ES I 

p scale diurnal variation all II 

p scale directional variation CC II 

p scale directional variation ES II 

p scale diurnal variation all III 

p scale directional variation ES III 

X shift all I 

y shift all I 

z shift all I 5 

X shift all II 

y shift all II 

z shift all II 5 

X shift all III 

y shift all III 

z shift all III 5 

X resolution all I 

y resolution all I 

z resolution all I 

X resolution all II 

y resolution all II 

z resolution all II 

X, y resolution constant term all III 0.065 

X, y resolution linear term all III —5.5 X 10""" 

X, y resolution quadratic term all III 3.9 X 10~^ 

z resolution constant term all III 0.0710 

z resolution linear term all III 1.16 X 10~ 

p resolution diurnal variation all I 

p resolution diurnal variation all II 

p resolution directional variation CC II 

p resolution directional variation ES II 



7 c- 



N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


+ 0.DD27 




scanned 


-0.0051 




In / rt. 


shirt 


and-refit 


N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


N/A 


shift 


and-reiit 


N /A 


shift 


and— refit 


N/A 




N/A 


N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


N/A 




N/A 


N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


1\ / rt. 


shift 


and-refit 


N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


N/A 


shift 


and-refit 


N/A 




N/A 



^ Correlation of -1. 
^ Correlation of -1. 

Correlation of ^ 

Correlation of 0.1 
Correlation of -1. 



1 -0.13 -0.74 \ 
0.13 1 0.31 
0.74 0.31 1 / 



Background neutrons from photo-disintegration observed 
with the PMT array in Phases I and II were weighted by 



(1 



,PMT 



+ en,corr.)(l + ^Pd)- 



(B22) 



This uncertainty was aheady included in the rates of 
these backgrounds in Table pCIX] for Phase III. 



Appendix C: Neutrino sensitivity 



Table XXII gives S{E^), the predicted spectrum of Ei, 
detectable by the SNO detector after including all effects 
such as the energy dependence of the cross-sections, re- 
action thresholds, and analysis cuts, but not including 
neutrino oscillations. 
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TABLE XVL cos^© PDF systematic uncertainties. 



Parameter Description 



Events Phase Nominal Variation 



Fit Value Application 





a„ 


COS 


50 


scale 


ES 


I 





±0.11 


N/A 


shift-and-refit 


e 
"o 


COS 


5© 


scale 


ES 


II 





±0.11 


N/A 


shift-and-refit 


e 
°-o 


COS 


50 


scale 


ES 


III 





±0.12 




scanned 


a" 


COS 


50 


scale directional variation 


ES 


I 





±0.022 


N/A 


shift-and-refit 


a" 

^dir 


COS 


50 


scale directional variation 


ES 


II 





±0.052 


N/A 


shift-and-refit 


a" 


COS 


50 


scale directional variation 


ES 


III 





±0.069 


-0 r)'\r. + 0.073 


scanned 



TABLE XVII. i?NCD PDF systematic uncertainties. 



Parameter Description 


Events 


Phase 


Nominal 


Variation 


^^NCD ^^^^ g^ai^, 

^^NCD E'ncd resolution 


n 
n 


III 
III 






±0.01 

+0.01 
-0.00 



TABLE XVIIL Background constraints in Phases I and IL The constraints were all applied to the combined day + night value. 



Background 



Phase Constraint 



126+*^ 

^ 1+1-4 

6.50 ± 1.11 

iQO+° °<'^ 

J^»u_0,054 



Fit Valu 
Day 



Night 



Application 



6.9 
+ 0.35 
0.34 
+3.3 
3.4 
+ 0-157 
0.157 
+ 0-059 
0.05^ 



[mBq] 
Internal ^"^Tl [mBq] 
External ^^^Bi [Bq] 
External ^'"^Tl [Bq] 
PMT [Arb.] 

AV surface neutrons [Arb.] 

AV =^"Bi [Arb.] 

AV ^"'^Tl [Arb.] 

hep neutrino [events] 

Other n [events] 

Atmospheric u [events] 

AV instrumental background [events] 

Internal '^^^Bi [Arb.] 

Internal ^"^Tl [mBq] 

Internal ^*'Na [mBq] 

External ^"Bi [Bq] 

External ^'"'Tl [Bq] 

PMT [Arb.] 

AV surface neutrons [Arb.] 

AV ^"Bi [Arb.] 

AV ^°**T1 [Arb.] 

hep neutrino [events] 

Other n [events] 

Atmospheric u [events] 

AV instrumental background [events] 



64.9 
1.11 
11.9 



'-7-1 
I +0.37 
-0-36 
3+4.2 
^-4-2 
n -I cQ + 0.202 
"■1^-^-0.199 
n Q'^S + 0-072 

0.9JS„o 
3.026tl;t?? 

Q ^00 + 2. 252 
z.ozz_2 

6 1 96+1-^1^ 
"■l»u_1.315 

N/A 
N/A 
N/A 
N/A 

"■ '^^-0.074 

69+^ ''*' 
274+'''3*2 

^■-^1^-0.342 

4.56 



+ 3.31 
3.35 
+ 0.159 



96.1 
1.09 
2.9 
0.265 
1.018 



n-216_o.i6o 
1 093+° 

n S91 +i-*S6 
fi 91O + 0.9S1 

N/A 
N/A 
N/A 
N/A 



N/r 
N/A 
N/A 
N/A 

-0.067 
-0.067 
,+ 1.41 
1-41 
-0.284 
-0.285 
+ 2.83 
■2-86 
U.UI i_o 133 

1 244+''-'''"' 



0.4957 
1.49! 

0.1931 
5.15'' 



N/5r 

N/A 
N/A 
N/A 



floated 
floated 
floated 
floated 
floated 
floated 
floated 
floated 
fixed 
shift-and-reflt 
shift-and-reflt 
shift-and-reflt 
floated 
floated 
floated 
floated 
floated 
floated 
floated 
floated 
floated 
fixed 
shift-and-reflt 
shift-and-reflt 
shift-and-reflt 



^ The fit was performed with day+nigiit combined, so there is only one fit value for both. 
^ Fixed at CC=0.35 SSM, ES=0.47 SSM, NC=1.0 SSM [21] . 
Correlation 1. 

One-sided, effect is symmetrized. 
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TABLE XIX. Background constraints in Phase III. The constraints were applied on the number of events observed in the NCD 
array. The number of events observed in the PMT array were obtained from the number of events in the NCD array multiplied 
by the PMT array ratio. 



Background 

External + AV 
Internal 
NCD bui; 
K2 hotspo" 
K5 liotspot 
NCD array cable: 
Atmospheric v aif 



NCD Array 
Constraint Fit Value 
40.9 ± 20.6 38.1 ± 19.2 



cosmogcnic muons 



31.0 ± 4.7 
27.6 ± 11.0 
32.8 ± 5.3 

8.0 ± 5.2 
13.6 ± 2.7 



30.9 ± 4.8 
27.2 ± 9.4 
32.6 ± 5.2 
45.4t',-3 

13.4 ± 2.7 



PMT Array 
ratio 
0.5037 
0.2677 
0.1667 
0.2854 
0.2650 
0.1407 
1.8134 



D /N asymmetry 



Constraint 
-0.020±0.011 
-0.034±0.112 








Fit Value 
-0.019±0.011 
-0.034±0.112 
N/A 
N/A 
N/A 
N/A 
N/A 



In the previous analysis of data from Phase III these two backgrounds were combined. 



TABLE XX. PMT background PDF parameters for the analytical PDF given by Equation [TO] in Section [IlL3l 



Parameter Phase Day /Night Constraint Fit Value Application 



e 




I day 




-6.73 ± 1.29 


-6 


„„+0.a5 

■^'^-0.56 




scanned 


e 




I night 




-5.64 ± 1.02 


-6 


40+0-31 




scanned 


e 




II day 




-6.26 ± 0.91 


-6 


70 + 0.29 
'O-0.37 




scanned 


e 




II night 




-6.98 ± 0.91 


-6 


79 + 0.24 
'^-0.33 




scanned 


f?] 


a 


I day 




Oil 


-0 


Y4+1.10 
'^-0.54 




scanned 


J?: 


5 


I night 




Oil 


-0 


QQ+0.39 
'-'^-0.12 




scanned 


f?: 


r 


II day 




Oil 





'^-0.26 




scanned 


f?: 


5 


II night 




Oil 





+0.26 
■^1-0.13 




scanned 


■n: 


i 


I day 




Oil 





nQ+0-62 
^^-0.61 




scanned 


m 


5 


I night 




Oil 





ns+0-77 




scanned 


v. 




II day 




Oil 


-2 


42+0.91 

^^-0.39 




scanned 


v. 


3 


II night 




Oil 


-3 


-., + 0.47 
'■^-0.49 




scanned 


Uo 


I day 





533 i 0.014 


n Kori +0.UU90 




scanned 


Wo 


I night 





533 i 0.014 


5469+° 




scanned 


Uo 


II day 





511 i 0.007 


n f;nnc + 0. 0055 
0.5096_o 0047 




scanned 


Uo 


II night 





511 i 0.007 


5119+" ™*^ 

U.011»_o 0055 




scanned 


U!l 


I day 





237 i 0.051 




N/A 


shift-and-rcfit 


OJl 


I night 





237 i 0.051 




N/A shift 


and-rcfit 


Ut 


II day 





182 i 0.095 




N/A 


shift-and-refit 


Ul 


II night 





182 i 0.095 




N/A 


shift-and-rcfit 


0. 


I day 





182 i 0.011 




N/A shift 


and-refit 


0. 


I night 





182 i 0.011 




N/A 


shift-and-refit 


p. 




II day 





195 i 0.007 




N/A 


shift-and-rcfit 


p. 




II night 





195 i 0.007 




N/A 


shift-and-refit 



b= -1.00 -I- 1.297?i, u = 6.63 + 0.93 * (0.60 * /?i + VT^^0Ji(fl-n2) 
^ b = 3.27+ 12.04r?i, u = 6.78-1- 1.52 * (0.96 * »yi -I- Vl - 0.96'^^) 
<: b = -0.33 -I- 2.08r7i, i/ = 5.32 + 1.01 * (0.91 * m + Vl - OM'^m) 

b = 0.49-f 3.02»yi, ly = 5.66 + 1.07* (0.94 * r;i + Vl - 0.942,72) 



TABLE XXI. Uncertainties in the neutron detection and photo-disintegration backgrounds. 



Parameter 


Description 




Phase 


Constraint 


Fit Value 


Application 


gPMT 


Neutron detection in the PMT 


array 


I 


± 0.0187 


N/A 


shift-and-rcfit 


gPMT 


Neutron detection in the PMT 


array 


II 


± 0.0124 


N/A 


shift-and-rcfit 


gPMT 


Correlated neutron detection in the PMT array 


I, II 


± 0.007 


N/A 


shift-and-rcfit 


gPMT 


Neutron detection in the PMT 


array 


III 


± 0.028 


-0.003 ± 0.028 


fioat 


gNCD 


Neutron detection in the NCD 


array 


III 


± 0.024 


-0.001 ± 0.023 


float 


epD 


Photo-disintegration 




I: II 


± 0.02 


N/A 


shift-and-refit 
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TABLE XXII. Predicted spectrum of E^, detectable by the SNO detector after including all effects such as the energy dependence 
of the cross-sections, reaction thresholds, and analysis cuts, but not including neutrino oscillations, S{Ev)- The number of 
events are for all three phases combined and assumes the BS05(OP) solar neutrino model ($b=5.69 x lO'^cm-^s"^) [21]. 



range CC ESe ES^^ range CC ESe ES^.^ 

[MeV] day night day night day night [MeV] day night day night day night 



2, 


.9- 


-3, 


.0 


0, 


.0 


0, 


.0 


0, 


,0 


0, 


.1 


0, 


.0 





.0 


9.1-9.2 




185, 


.3 


231 


.3 


19, 


.9 


24 


.8 


3 


.2 


4, 


.0 


3, 


.0- 


-3, 


.1 


0, 


.0 





.0 


0, 


,1 


0, 


.1 


0, 


.0 





.0 


9.2-9.3 




187, 


.7 


235 


.1 


20, 


.1 


24, 


.6 


3 


.2 


4, 


.0 


3, 


.1- 


-3, 


.2 


0, 


.0 





.0 


0, 


,1 


0, 


.1 


0, 


.0 





.0 


9.3-9.4 




191 


.3 


238 


.8 


20, 


.0 


24 


.7 


3 


.2 


4, 


,0 


3, 


.2- 


-3, 


.3 


0, 


.0 





.0 


0, 


,1 


0, 


.2 


0, 


.0 





.0 


9.4-9.5 




194, 


.0 


241 


.2 


20, 


.1 


24 


.8 


3 


.2 


3, 


,9 


3, 


.3- 


-3, 


.4 


0, 


.0 





.0 


0, 


,2 


0, 


.2 


0, 


.0 





.0 


9.5-9.6 




196, 


.4 


244 


.4 


19, 


.7 


24, 


.4 


3 


.2 


3, 


.9 


3, 


.4- 


-3, 


.5 


0, 


.0 





.0 


0, 


,3 


0, 


.3 


0, 


.0 





.1 


9.6-9.7 




198, 


.6 


246 


.1 


19, 


.9 


24 


.5 


3 


.2 


3, 


,9 


3, 


.5- 


-3, 


.6 


0, 


.0 





.0 


0, 


,4 


0, 


.5 


0, 


.1 





.1 


9.7-9.8 




199, 


.7 


249 


.3 


19, 


.5 


24, 


.2 


3 


.1 


3, 


.9 


3, 


.6- 


-3, 


.7 


0, 


.0 





.0 


0, 


.5 


0, 


.6 


0, 


.1 





.1 


9.8-9.9 




199, 


.8 


250 


.2 


19, 


.4 


24, 


.0 


3 


.1 


3, 


.9 


3, 


.7- 


-3, 


.8 


0, 


.0 





.0 


0, 


,6 


0, 


.7 


0, 


.1 





.1 


9.9-10.0 


200, 


.1 


251 


.1 


19, 


.3 


23, 


.7 


3 


.1 


3, 


,8 


3, 


.8- 


-3, 


.9 


0, 


.0 





.0 


0, 


,8 


0, 


.9 


0, 


.1 





.2 


10.0-10, 


.1 


200, 


.9 


251 


.0 


19, 


.0 


23, 


.5 


3 


.0 


3, 


.7 


3, 


.9- 


-4, 


.0 


0, 


.0 





.1 


0, 


,9 


1 


.2 


0, 


.2 





.2 


10.1-10, 


.2 


200, 


.8 


250 


.8 


18, 


.7 


23, 


.1 


3 


.0 


3, 


,7 


4, 


.0- 


-4, 


.1 


0, 


.1 


0, 


.1 


1 


,1 


1 


.5 


0, 


.2 





.2 


10.2-10, 


.3 


199, 


.7 


248 


.7 


18, 


.5 


22, 


.9 


2 


.9 


3, 


.6 


4, 


.1- 


-4, 


.2 


0, 


.1 





.2 


1 


,4 


1 


.6 


0, 


.2 





.3 


10.3-10, 


.4 


198, 


.3 


246 


.8 


18, 


.1 


22, 


.3 


2 


.9 


3, 


.6 


4, 


.2- 


■4, 


.3 


0, 


.3 





.3 


1 


,6 


2, 


.0 


0, 


.3 





.3 


10.4-10, 


.5 


198, 


.0 


246 


.5 


17, 


.7 


21 


.9 


2 


.9 


3, 


,5 


4, 


.3- 


-4, 


.4 


0, 


.4 





.5 


1 


,9 


2, 


.4 


0, 


.3 





.4 


10.5-10, 


.6 


195, 


.4 


244 


.5 


17, 


.3 


21 


.4 


2 


.8 


3, 


.4 


4, 


.4- 


-4, 


.5 


0, 


.6 





.8 


2, 


,2 


2, 


.7 


0, 


.4 





.4 


10.6-10, 


.7 


193, 


.9 


241 


.6 


17, 


.0 


21 


.1 


2 


.7 


3, 


,4 


4, 


.5- 


-4, 


.6 


1 


.0 


1 


.2 


2, 


.5 


3, 


.1 


0, 


.4 





.5 


10.7-10, 


.8 


191 


.7 


237 


.5 


16, 


.5 


20, 


.6 


2 


.6 


3, 


.3 


4, 


.6- 


-4, 


.7 


1 


.4 


1 


.8 


2, 


,8 


3, 


.5 


0, 


.5 





.6 


10.8-10, 


.9 


187, 


.1 


233 


.2 


16, 


.2 


19, 


.8 


2 


.6 


3, 


.2 


4, 


.7- 


-4, 


.8 


2, 


.0 


2 


.7 


3, 


,1 


3, 


.9 


0, 


.5 





.6 


10.9-11 


.0 


184, 


.4 


229 


.6 


15, 


.9 


19, 


.4 


2 


.5 


3, 


,1 


4, 


.8- 


-4, 


.9 


2, 


.8 


3, 


.6 


3, 


,5 


4, 


.5 


0, 


.6 





.7 


11.0-11 


.1 


180, 


.5 


226 


.4 


15, 


.3 


18, 


.8 


2 


.4 


3, 


.0 


4, 


.9- 


-5, 


.0 


3, 


.8 


4, 


.8 


3, 


,9 


4, 


.8 


0, 


.6 





.8 


11.1-11 


.2 


176, 


.6 


220 


.0 


14, 


.7 


18, 


.3 


2 


.3 


2, 


,9 


5, 


.0- 


-5, 


.1 


4, 


.9 


6, 


.3 


4, 


,4 


5, 


.3 


0, 


.7 





.9 


11.2-11 


.3 


172, 


.2 


214 


.3 


14, 


.2 


17, 


.5 


2 


.3 


2, 


.8 


5, 


.1- 


-5, 


.2 


6, 


.4 


8, 


.0 


4, 


,7 


5, 


.9 


0, 


.8 


1 


.0 


11.3-11 


.4 


167, 


.8 


208 


.2 


13, 


.6 


16, 


.8 


2 


.2 


2, 


.7 


5, 


.2- 


-5, 


.3 


8, 


.0 


10, 


.2 


5, 


,2 


6, 


.4 


0, 


.8 


1 


.0 


11.4-11 


.5 


162, 


.4 


202 


.7 


13, 


.1 


16, 


.3 


2 


.1 


2, 


,6 


5, 


.3- 


-5, 


.4 


10, 


.0 


12, 


.7 


5, 


,7 


6, 


.9 


0, 


.9 


1 


.1 


11.5-11 


.6 


157, 


.5 


195 


.9 


12, 


.6 


15, 


.6 


2 


.0 


2, 


.5 


5, 


.4- 


-5, 


.5 


12, 


.3 


15, 


.5 


6, 


,1 


7, 


.6 


1 


.0 


1 


.2 


11.6-11 


.7 


152, 


.2 


188 


.9 


12, 


.0 


14, 


.9 


1 


.9 


2, 


,4 


5, 


.5- 


-5, 


.6 


14, 


.5 


18, 


.7 


6, 


,6 


8, 


.1 


1 


.0 


1 


.3 


11.7-11 


.8 


146, 


.5 


181 


.8 


11 


.4 


14, 


.2 


1 


.8 


2, 


.3 


5, 


.6- 


-5, 


.7 


17, 


.4 


22, 


.1 


7, 


,1 


8, 


.8 


1 


.1 


1 


.4 


11.8-11 


.9 


139, 


.9 


173 


.9 


11 


.0 


13, 


.4 


1 


.7 


2, 


.1 
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